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Abstract 

We present an in-depth analysis of the concept of spin precession frequency for 
integrable orbital motion in storage rings. Spin motion on the periodic closed orbit of 
a storage ring can be analyzed in terms of the Floquet theorem for equations of motion 
with periodic parameters and a spin precession frequency emerges in a Floquet expo- 
nent as an additional frequency of the system. To define a spin precession frequency 
on nonperiodic synchro-betatron orbits we exploit the important concept of quasiperi- 
odicity. This allows a generalization of the Floquet theorem so that a spin precession 
frequency can be defined in this case too. This frequency appears in a Floquet-like 
exponent as an additional frequency in the system in analogy with the case of motion 
on the closed orbit. These circumstances lead naturally to the definition of the uni- 
form precession rate and a definition of spin tune. A spin tune is a uniform precession 
rate obtained when certain conditions are fulfilled. Having defined spin tune we define 
spin-orbit resonance on synchro-betatron orbits and examine its consequences. We 
give conditions for the existence of uniform precession rates and spin tunes (e.g. where 
small divisors are controlled by applying a Diophantine condition) and illustrate the 
various aspects of our description with several examples. The formalism also suggests 
the use of spectral analysis to "measure" spin tune during computer simulations of spin 
motion on synchro-betatron orbits. 
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1 Introduction 

This paper provides a rigorous discussion of the concept of spin precession frequency on 
synchro-betatron orbits in storage rings. To set the scene we begin by introducing some 
key physical ideas via the equations of orbit and spin motion and the notion of spin-orbit 
equilibrium. 

The spin expectation value S ("the spin") in the rest frame of, for example, a proton, an 
electron or a muon moving in electric and magnetic fields under the influence of the Lorentz 
force precesses according to the Thomas-Bargmann-Michel-Telegdi (T-BMT) equation pQ 

dS/dt = 0(E, B, v) x S , (1.1) 

where the precession vector ft depends on E, B and v which are respectively the electric and 
magnetic fields, and the velocity. Particle motion with respect to the synchronous closed, i.e. 
periodic, orbit is described in terms of three pairs of canonical variables which we combine 
into a vector u with six components. For example, two of the pairs can describe transverse 
motion and one pair can describe longitudinal (synchrotron) motion within a bunch. Since 
we are dealing with storage rings we take the orbital motion to be bounded. In this paper 
we ignore radiation, interparticle interactions and interactions with the vacuum system. In 
(11.1)1 the independent variable is the time t. However, since the electric and magnetic guide 
fields in particle accelerators and storage rings are fixed in space it is more convenient to 
adopt the standard practice of replacing t with the angular position around the ring, the 
azimuth 9 = 2tts/L, where s is the distance around the ring and L is the circumference. 
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Since E, B and v depend on u and 9 we can now rewrite (jl.lj) in the form 

dS/d9 = ft(9,u(9j)xS (1.2) 

where f2 is the precession vector obtained from ft by rescaling with dt/d9 and transforming 
to machine coordinates 0. 

If the beam is in equilibrium, i.e. if the the phase space density p(9, u) is 27r-periodic 
in 9, we then write it as p eq so that p eq (9,u) = p eq (9 + 2n,u). We normalize it to unity: 
J dup eq = 1. The necessary condition for this kind of equilibrium, namely that at a fixed u 
the fields are 27r-periodic in 9, is automatically fulfilled in a storage ring. But of course the 
boundedness of the motion is also required. Conditions necessary for beam equilibrium and 
a way of calculating p eq using ergodic theory and the concept of "stroboscopic averaging" 
are described in detail in 

The statistical properties of the spins are encoded in the quantum mechanical spin density 
matrix. But for spin 1/2 particles this can be completely parametrized by the polarization 
vector For particle beams we need the local polarization P\ oc (9,u) at each point in 
phase space u. P\ oc {9,u) is the average of the normalized spin vectors, S/|S|, at u, where 
| • | denotes the Euclidean norm. We define the polarization of the whole beam, the "beam 
polarization", at a given azimuth as f dup eq P\ oc . 

Since the T-BMT equation (jl.2j) is linear in S and since the particles at (9, u) all see the 
same Q(9,u), P\ oc (9,u(9)) also obeys the T-BMT equation Furthermore, the length 
of P loc (9, u) is constant along a phase space trajectory. For a storage ring at fixed energy, 
O is 27r-periodic in 9 at a fixed position in phase space u so that fl(9,u) = £1(9 + 27t,m). 
This opens up the possibility of a spin distribution that is the same from turn to turn, i.e. 
in equilibrium. Then P\ oc (9,u(9)) not only obeys the T-BMT equation, but P\ oc (9, u) is 
27r-periodic in 9 for fixed u and we then write it as P eq so that P eq (6',M) = P eq (6' + 2tt,u). 
We denote the unit vector along P cq (9, u) by n(0, u). This also obeys the T-BMT equation 
along orbits and is 27r-periodic in 9: n(9, u) = n(9 + 2n, u). 

The method described in |3] for constructing p eq can be extended as in jB] for constructing 
n and from the treatments in j3] it is clear that the existence of p eq and n do not require 
that the orbital motion be integrable. But, of course, the conclusions of [3] are still valid if 
the motion is integrable. Moreover, particle motion in storage rings is usually close enough 
to integrability to allow the motion to be characterized in terms of well defined betatron and 
synchrotron frequencies. This, in turn, allows predictions to be made about beam stability 
via the concept of orbital resonance. Thus, in the remainder of this paper we will assume 
that the orbital motion is integrable. Then, as we shall see, the stability of spin motion can 
also be discussed in terms of resonance, namely "spin-orbit resonance". Of course, integrable 
orbital motion and spin-orbit equilibrium are idealizations. Nevertheless, these idealizations 
often provide useful starting points for calculations. 

For integrable particle motion the position of a particle in phase space is represented by 
three pairs of action-angle variables ( Jj, <pi, i = 1,2, 3) and is determined by a Hamiltonian 
H(J). Thus the orbital phase space is partitioned into disjoint tori, each of which is char- 
acterized by a unique set of Js. We now define u := (0i, 02, 4>s, Ji, J2, J3) = (0, J)- The 
actions are constants of the motion and for fixed J the constant rate of advance of each 
LUi(J) := dH/dJi = d(pi/d9, is called an orbital tune. These frequencies are the number of 
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oscillations per turn around the ring. In beam physics such frequencies are often referred to 
as tunes and we have adopted that usage. We will only consider storage rings running at 
fixed nominal energy. 

For integrable motion, the 27r-periodicity in 9 of Q, P eq and n is accompanied by 2n- 
periodicity in 0i,02 and 03. So as well as being a solution to the T-BMT equation along 
orbits u(6), n satisfies nontrivial periodicity conditions. In our later discussions on quasiperi- 
odicity we will require that it also depend sufficiently regularly on the azimuth and the orbital 
angles, i.e. n must be "smooth" in the sense defined in the main text. This corresponds to 
the expectation that P eq also be smooth. The equilibrium density p eq is also 27r-periodic in 
9 and off orbital resonance it just depends on J. 

Since for every J the field n(9, 0, J) is invariant from turn to turn, it is now often called the 
invariant spin field (ISF). The ISF is a central object in the theory of polarization in storage 
rings [II IHJ IIE E3 • F° r example, for an ISF sufficiently regular in J, off orbital resonance and 
away from the spin-orbit resonances to be defined below, an upper limit to the equilibrium 
beam polarization at a particular 9 is J dJp eq (J)\ J d<p n(6, 0, J)| and it is reached only when 
the J dcf) n are parallel. This is easy to see by noting that if |P eq (#, 0, J)| were to vary over 
a torus, the beam polarization would vary from turn to turn. So equilibrium implies that 
|P eq (#, 0, J) | is constant over a torus. The maximum equilibrium polarization on each torus 
is reached when |P eq (6 l , 4>,J)\ = 1- Note that a zero value for J dJp eq (J)\ J dtp n| at some 9 
does not mean that the beam is depolarized. It could well be that the beam is fully polar- 
ized at each point in phase space but that the geometry of n causes the integral to vanish. 
Then, if a change of parameters were to change the geometry of n so that the integral were 
to become nonzero, and if the change were carried out adiabatically, the beam polarization 
would reappear. Furthermore, the fact that the integral vanishes at one position 9 does not 
mean that it vanishes at other positions. We prefer to reserve the term "depolarization" for 
a definitely irreversible loss of polarization such as occurs in the presence of noise or for an 
effectively irreversible loss of the kind that can occur when spin-orbit resonances are crossed 
See Sec. X too. Although we have introduced n via the notion of equilibrium, the inte- 
gral J dJp eq (J) \ J d<j) n| also contains useful information when the spin distribution is not in 
equilibrium: J dJp eq (J)\ J d(fi n| gives an upper limit for the time averaged polarization away 
from spin-orbit resonances. The maximum is reached on each torus when the polarization 
is in equilibrium and with \P eq (9, 4>, J)\ = 1- See O Section 2.2.8] and [T3J Section 4.4]. An 
example of the origin and behavior of nonequilibrium beam polarization is given in Figure 
9 in [B] where large oscillations are evident. However, polarimeters and particle detectors 
cannot collect data quickly enough to make such oscillations observable. Instead, only the 
time averaged polarization can be observed or exploited. But as we have just seen, we can 
still estimate its maximum value. That depends only on the geometry of n and it is reached 
for each torus when the spread of n is minimized. The ISF also provides a perfect tool for 
estimating the long term effects on the beam polarization of small perturbations such as 
radiation [7] or electric and magnetic fields which cause nonintegrable orbital motion. In 
particular, one begins with a spin-orbit system which is invariant from turn to turn, i.e. with 
an equilibrium orbital distribution and with spins set parallel to the n(9,u). Then, since 
the system is initially in equilibrium, the effects of the perturbations cannot be masked by 
the natural, potentially large, variations of the beam polarization of the kind depicted in 
Figure 9 in 0. An ability to construct n(9,u) for integrable orbital motion and understand 
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its behavior is then indispensable. 

For our integrable orbital motion of electrons and protons and up to energies of a few 
GeV, an approximate n can be calculated in a first order perturbation theory by an extension 
of the code SLIM [T^ITH]. and in higher order perturbation theory by the codes, SMILE 
Forget-me-Not ^B] and SpinLie However, for the high magnetic fields characteristic of 
proton rings running at energies of hundreds of GeV, perturbative methods are inadequate. 
Then the method of stroboscopic averaging as in the code SPRINT jH] should be used. This 
is a numerical, nonperturbative algorithm and yields high accuracy for real rings even when 
all modes of orbit oscillation are included simultaneously. One can also use Fourier methods 
as in the codes SODOM2 03j or MILES [IH]. SODOM2 has been very useful for orbital 
motion restricted to one plane [13J. MILES gives explicit formulae which are applicable to 
some simple models. So far, the only practical general way to calculate the invariant spin 
field is to use stroboscopic averaging. 

As for any dynamical system we hope to understand more about spin motion by studying 
its spectrum of frequencies. Various quantities, which seem at first sight to be related to 
spin frequencies, can be found in the literature and we will mention some in Section 10. But 
a true component of a spectrum quantifies long term behavior. Thus any definition of spin 
precession frequency should reflect that stipulation. The choice can be further narrowed 
by requiring that the spectrum give useful clues about the behavior of sets of spins, and 
in particular about the beam polarization. After all, the experimenters using the beams in 
storage rings are just interested in the beam polarization, not individual spins. 

Experience has shown that the best choice for characterizing spin motion in storage rings 
is the traditional one |2"U1 |2~T| |2"2"] . namely the so-called amplitude dependent spin tune (briefly 
"spin tune"), which we usually denote by uj s . Assuming n exists, the spin tune measures the 
number of spin precessions around n(9,u), per turn around the ring, for a particle on the 
orbit u{9) and it provides a way to quantify the correlation between the spin motion and the 
orbital motion which "drives" it, and thereby forecast a qualitative aspect of spin motion, 
namely the degree of regularity of the spin motion. In particular, the spin motion can in 
general become very erratic when a spin tune is near a low order resonance condition 

uj s = m + m\ LUi + m 2 iO% + m 3 u 3 = m ■ (1, lu), (1.3) 

where m is a vector of integers and the quantity \mx\ + \m 2 \ + |m 3 | is usually called the 
order of the resonance. Correspondingly, close to spin-orbit resonance n(6, u) can become 
a very sensitive function of u. This sensitivity has immediate consequences for work with 
polarized beams. For example, the maximum attainable equilibrium beam polarization of 
a stored high energy proton beam can be unacceptably low or the rate of depolarization, 
due to synchrotron radiation, of a stored electron beam can be unacceptably high [7] . Note, 
however, that the \fd(j) n|/(27r) 3 on a torus can sometimes be small away from spin-orbit 
resonance and that proximity to a spin-orbit resonance, especially one of high order, does not 
automatically imply that the \J d(f) n|/(27r) 3 on a torus is low. The resonance might be very 
weak. Another feature of our definition of spin frequency is that, as we shall see, it is this 
quantity whose spectrum one obtains in a straightforward spectral analysis of spin motion 
during spin-orbit tracking simulations. In other words: in an ideal world with technology 
which could select particles on a torus at a fixed J, it could be measured. 
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Right at the resonance condition (jl.3|) . n is in general nonunique. However, as we shall 
see, our spin-orbit systems exhibit a tendency to avoid exact spin-orbit resonance. Since 
Q in general depends on J and the particle energy on the closed orbit, u s usually varies 
with J (hence "amplitude dependent spin tune"), and the particle energy on the closed 
orbit. We will call the latter the "beam energy". We emphasize that n is a field over 
the six dimensional phase space so that synchrotron motion is built in from the start. Thus 
although uj s varies with the beam energy and Ji, J 2 and J3, it does not change during a period 
of synchrotron motion. If n were defined on four dimensional transverse phase space and 
the energy oscillations due to synchrotron motion were added as an afterthought, it would 
not be useful for describing equilibrium polarization. Instead, we would have to characterize 
the beam polarization using time averages. We return to this theme in Section 10. On the 
closed orbit, i.e. for J x = J 2 = J3 = 0, an n exists which is independent of 0. We denote it 
by n 

The calculation of spin tune on the closed orbit presents no problem: it can, as we shall 
see, be extracted from an eigenvalue of the 1— turn spin map. But the definition of spin 
tune for J ^ 0, i.e. on synchro-betatron orbits is much more subtle. Moreover, it requires 
precision. Notions of spin frequency for synchro-betatron orbits appearing in the literature 
are often not precisely presented and some appear to possess no capacity for predicting the 
qualitative aspects of spin motion in storage rings. 

This brings us to the purpose of this paper. This is to provide a rigorous discussion of 
the concept of spin precession frequency on synchro-betatron orbits and thereby consolidate 
a framework for systematizing and classifying spin motion in storage rings. For this we make 
a careful mathematical study of the consequences for spin motion of the periodicities of ft 
in 9,4>i,4>2 and 03, using precise definitions and carefully formulated theorems and we make 
use of the ISF and other concepts which we distill from the literature and "folklore" on spin 
dynamics in storage rings [SI [3 El El H21 CHI CHI HH1 1201 IZH [221 [23 I2H [213 For example, 
we will show that under the appropriate conditions, the existence of the ISF with the above 
mentioned periodicities implies that the 9 dependence of the general solutions of (jl.2|) will 
contain five frequencies. Four of them are the orbital tunes Ui,u 2 ,u) 3 and the circulation 
tune u c — 1, i.e. the frequency associated with the 27r-periodicity in 9. A fifth tune emerges 
which, under circumstances to be described, is a spin tune uj s . The general solutions will then 
be found to be quasiperiodic with the tunes 1, u>i, cu 2 , u 3 , uj s . Moreover the results obtained 
here can be viewed as a generalization of Floquet theory. Given the confusion surrounding 
definitions of spin precession frequencies, the treatment of the kind that we provide here 
seems to be very necessary. Our assumptions about O are weak enough to cover the situations 
of most interest for storage rings, namely typical integrable synchro-betatron motion. Several 
of our theorems assume the existence of n but although we have ways to find approximate 
n, the determination of complete conditions for its existence is an outstanding mathematical 
issue. This question can, for example, be investigated using ergodic theory [23 [2Zj and the 
method of stroboscopic averaging. See [3]. Moreover, simulations indicate that approximate 
ISFs do exist. This means that one obtains objects which, at least approximately, behave 
like an ISF. Moreover in some instances approximations even lead to an Q, for which an exact 
ISF can be found, e.g. as in the single resonance model - see Section 7. 

Although we have introduced the vector n by studying spin-orbit equilibrium, it was 
first discussed by Derbenev and Kondratenko as a vehicle for constructing joint action- 
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angle variables for spin and orbital motion from their semiclassical spin-orbit Hamiltonian 
[2~1"| 123 121 HI E3- This Hamiltonian is derived from the Dirac Hamiltonian by a Foldy- 
Wouthuysen transformation taken to first order in h jZj. In that picture the spin tune 
emerges as the rate of advance of a spin phase [2H 1221 12] • The terms at first order in h in the 
Derbenev-Kondratenko Hamiltonian are those containing spin and these terms imply a force 
of the Stern-Gerlach (S-G) type 4J. The S-G forces on trajectories appear at first order 
in h and a "back reaction" on the spin of the S-G perturbation to the orbit would involve 
an addition to the spin precession rate of order h. However, in this paper the effect of S-G 
forces on spin and orbit motion is neglected and we just operate with the Lorentz force and 
the T-BMT equation and for the Lorentz force and ft just include, as is usual, the terms of 
zeroth order in ft. There are several reasons for this approach. First, it is far from clear what 
form the S-G forces should take. In fact there is considerable ambiguity in the choice of the 
S-G forces. This is covered in detail in [2E]- See "2H] too and the bibliography in [211 • The 
second ground has to do with the size of the S-G forces. Since the S-G forces are of first 
order in h they are extremely small compared to the Lorentz forces which are of zeroth order 
in H . They are also small compared to typical spurious perturbations to trajectories like 
noise and collective effects. So S-G forces would not cause changes of practical significance 
to the results that we present. In particular, in practical situations in a storage ring there 
would be no significant change in the phenomenology of spin-orbit resonances even if the 
S-G forces were to cause tiny changes in the orbital tunes. The third ground is that it is 
far from clear that it makes sense to treat an essentially quantum mechanical system with 
a classical "over interpretation" of the influence of the S-G force on the spin. An example 
of an effect which is not taken into account by a naive application of classical S-G forces, is 
given in [30 . It is implied there that long term shifts of an orbit due to S-G forces will be 
nullified when the spin undergoes a quantum flip and the S-G force then acts in the reverse 
direction. See [21 p. 137] for a classical perspective on this. In summary, we believe that a 
too literal interpretation of the S-G-like forces in the semiclassical spin-orbit Hamiltonian 
could lead to manipulations and conclusions of little relevance and utility for illuminating 
the core phenomenology of spin motion in typical storage rings [3T1 13*2] . We believe that 
the first priority is to begin with just the Lorentz force and the T-BMT equation. Then, as 
mentioned earlier, once the equilibrium state of the system has been defined, other influences 
such as nonlinear fields, noise, collective effects, synchrotron radiation and the very small 
S-G-like effects can be included as perturbations. 

The paper is structured as follows. In Section 2 we begin by discussing some important 
consequences of (|1.2j) . Here we introduce the central concept of a uniform precession frame 
(UPF) and the associated uniform precession rate (UPR). The UPF provides a coordinate 
system for spin. Then in Section 3 we give a detailed discussion of spin motion on the 
closed orbit where ft is independent of <fr so that the Floquet theorem applies. Sections 2 
and 3 contain standard results but we present them in forms which motivate their extension 
in later sections. Section 4 contains the definition of a quasiperiodic function and collects 

1 For example, at the r.m.s. radius of the 920 GeV/c proton beam in a quadrupole magnet in the 
HERA proton ring |12l ITS] , a rough estimate for the transverse S-G force from the Derbenev-Kondratenko 
Hamiltonian gives a S-G force which is about 10~ 12 of the Lorentz force on a proton. At a fixed radius this 
ratio is essentially independent of the beam energy. The S-G energy at that radius is of the order of 10 -21 of 
the kinetic energy. The S-G energy in a HERA dipole magnet is of the order of 10 -19 of the kinetic energy. 
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some properties useful for the discussion following. In particular it defines a Diophantine 
condition needed for handling a problem with small divisors. The key ideas are formalized 
in Lemmas 4.3, 4.7 and 4.8. Section 5 uses the concept of a UPF, quasiperiodic with orbital 
frequencies, to define the proper UPR, the spin tune and spin-orbit resonance. The main 
theorem in Section 5 is Theorem 5.3 which allows us to define equivalence classes of spin 
tunes. The presentation in Sections 2, 3, 4 and 5 is deliberately rather general and abstract. 
Then in Section 6 we introduce a field called the invariant frame field (IFF) which is used to 
construct UPFs. There we consider the angular phase space as a whole to prove theorems 
about the concepts introduced in Section 5. We also connect the abstract ideas introduced 
earlier to a familiar physical idea, namely that if the orbital tune were off orbital resonance 
(m • (l,o;) = only with the vector of integers m — 0), the existence of a nonunique ISF 
would imply that the system were on spin-orbit resonance. The main theorems in Section 
6 are Theorems 6.3 - 5. Theorem 6.3 is used in the proof of Theorem 6.5 and it is generally 
our main tool for showing that a torus is "well-tuned". The proof of Theorem 6.3 relies on 
Theorem 5.3. Some examples of the formalism for model fts are presented in Sections 6, 7 
and 8. Note that, except for some examples, we allow the number of action-angle pairs, d, 
to be arbitrary (but > 1) although for spin motion in storage rings, the case d — 3 is the 
most important. To aid the reader we mark the key equations with a *** on the left. 

As a byproduct of the quasiperiodic structure of the solutions we suggest using spectral 
analysis as a way of "measuring" the spin tune during spin-orbit tracking simulations and 
thereby complementing other methods already in use |BJ EH HH1 EEZI ECU El EH]- Spectral 
analysis may also lead to a practical method for deciding whether an invariant spin field 
exists. These ideas are presented in Section 9 and formalized in Theorems 9.1 and 9.2. 

The paper is summarized in Section 10 where our concepts are also related to simulations 
and used to discuss some popular notions. 

For the rest of the paper, apart from Section 10, we will now adopt a more efficient 
notation whereby we use the symbols J = (J\, Jj), uj = (u\, ...,ujd) and = (0i, ...,0d), 
(J,w,<j> G R d ) to mean respectively the list of orbital actions, orbital tunes and orbital 
angles. From now on we will also adopt the frame dependent abbreviations S = (S±, S 2 , S 3 ), 
O = (Qi, 2 , 3 ) and n = (nj, n 2 , n 3 ). Generally, if appears as an independent variable in 
a function, the function will be 27r-periodic in <fid- I n that case we say for brevity that 

the function is 27r-periodic in 0. 

In terms of the new notation, the T-BMT equation and the equations of orbital motion 

are 



where A is a real skew-symmetric 3x3 matrix with nonzero elements Au = —A21 = 
—0 3 , A13 = —A31 = 2 and A23 = —A32 = —Oi. The dot over a symbol denotes differen- 
tiation w.r.t. 9. Because the J dependence is only parametric, we will often suppress the 
symbol J in u and A, e.g. as in A(9, 0). Clearly, A is 27r-periodic in 9 and in 0. For brevity 
we just say that functions with such periodicity are 27r-periodic. On the closed orbit, i.e. 
for J = 0, A is independent of 0. Note that on the torus J = the angular variables 
play a largely artificial role because here A is independent of 0. But their inclusion is very 



S 



J 



A(9,<p, J)S, 
0, = w(J), 



(1.4) 
(1.5) 
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convenient as it allows one to treat all tori on the same basis. Then all definitions, e.g. that 
of the ISF, apply to all tori. 

A function is called C r , if the function together with all of its partial derivatives up to 
and including those of order r are continuous. In this paper we will assume that, for fixed 
J, A is a C 1 function of (9,(f>). A C 1 function will be called smooth. The smoothness of 
A corresponds to the fact that in real storage rings, the magnetic and electric fields are 
smooth functions of space and time. The labels for the definitions, propositions, theorems 
and lemmas are chosen in a way which indicates their relative positions in the text. 



2 General Properties of the Spin Motion 

We begin by establishing some basic components of our formalism. 

Clearly (jl.5J) gives J{9) = J and <p(9) = u(J )9 + <p Q where J and O are the actions 
and phases at 9 = 0. Thus an orbit is labeled by (<fto, Jo). But if we consider a fixed torus 
Jo we often suppress the symbol Jo- By "a fixed torus Jo" we mean that the orbital tune 
has the value u( Jo) and that the spin motion is characterized by the function A(9, <f), Jq) of 
9 and (p. Equation (jl.4|) thus becomes 

**k S = A(9;(j) ,uj)S , S(0; fa, u) = S , (2.1) 

where the real skew-symmetric A is defined by 

A{9]fa,uo) :=A{9,u9 + <j) ) . (2.2) 

As will become clear from Definition 4.1 in Section 4, A{9] 4>o,uj) is a quasiperiodic function 
of 9 with the tunes (frequencies) 1, U\, Ud- The solutions of (j2.1j) can be written as 
S(9; 4>o,uj) = $(#; 4>o, co)S(0] (fio, to) in terms of the principal solution matrix at (f> which is 
the 3x3 matrix $(6 I ; O , u), defined uniquely by the initial value problem 

= A(0; M *, (2.3) 
4(0; fo,u) = /. (2.4) 

Thus the principal solution matrix at 0o is the spin transport matrix from the azimuth to 
the azimuth 9. Occasionally we call a solution S(9; O ) of (|2.1|) a "spin trajectory" at 4> Q . 
The choice 9 = for the starting azimuth does not imply a loss of generality as can be seen 
by considering the general initial value problem S^o) = So, <K#o) = 4>o for (jl.4j) and (jl.5j) . 

Note that A(9; <p ) and 0o) are 27r-periodic in O an d that by the smoothness of A 
the principal solution matrix is a smooth function of (9, (fro) |34^ I35j. We will sometimes 
suppress the symbols 4> and uj in S, A and $. The key property of $(#) is that it belongs 
to SO(3), i.e. 

$ T (#)$(#) = / = $(#)$ T (#) and det($(0)) = 1 , (2.5) 

as is easily proved using (|2.3j) and (|2.4|) . 

Let Sq and Sq be two initial conditions for (j2.1j) and let a ■ b = a\bi + 02^2 + 0363 be the 
real inner product. Then S\9) ■ S 2 (9) = (^(9)S^) ■ ($(0)S O 2 ) = ($ T (0)$(0)S' O 1 ) • S 2 = 5 x • S 2 
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so that the inner product of any two solutions of (|2.1|) is conserved. In particular, the length 
of a spin vector and the angle between any two spin trajectories at the same 0o is conserved. 
In addition, it is easy to show that the cross product of two solutions is a solution. In the 
remainder of this paper we will, for convenience, allow spins S to have arbitrary length. 

An interesting property of (|2.1j) is that knowledge of one solution completely determines 
$ by a simple integration. It is a standard result for linear systems that knowledge of one 
solution can be used to reduce the dimension by one. Here it reduces the dimension by two 
because of the special structure of $ in ()2.5|) as we will now demonstrate. 

Let v 3 {9) be a solution of (|2.1j) . i.e. a spin trajectory at 0o, and let v 3 be of norm 1. 
Choose v l {9) and v 2 (Q) so that 

V{0):= [v\9),v 2 (9),v 3 (9)} (2.6) 

is a 5*0(3) matrix. One can for example require v 1 and v 2 to be solutions of v k = (v 3 x v 3 ) x v h 
(k=l,2), whence we can assume that V(9) is smooth, i.e. a C 1 function. 
Next we make a transformation $ — > \1/ on (J2.3)) and (|2.4|) defined by 

$ = V(0)V. (2.7) 

This gives 

* = C(0)tf, tf(0) = ^ T (0) , (2.8) 

where 

C{9) = V T {9) (a{9)V{9) - V{6)} . (2.9) 

Since V{9) e S0{3), V T {9)V{9) = -V T {9)V{9) and thus C T = V T {A T V + V) = -C by the 
skew symmetry of A. Therefore C is skew-symmetric as expected for rotations. The third 
column of C is 

/ \ 

= V T (Av 3 - v 3 ) = 0, 

since v 3 is a solution of (|2.1j) . and the skew symmetry of C yields 

/ -1 

C(9) = c v (9) 1 | =:c v (9)J, (2.10) 
\ 

where fl2.10|l also serves to define J . Therefore 




tf(0) = exp \ J J c v (9')d9' | V 1 (0) 
as is easily checked by differentiation. Finally, 

$(0; (f> ) = V(9) exp J J c v (9')d9' J V T (0). (2.11) 
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The exponentials in these equations can be evaluated by noting that 

(cos r — sin r \ 
sinr cost . (2.12) 
1/ 

So we have constructed the complete principal solution matrix by starting from just one 
solution of ()2.1|) and assuming the existence of a smooth V. This is the result we were aiming 
for. 

Note that by flZHl) , (GUIl) we have 

V(9) = A(9; <j> )V(6) - c v (9)V(9)J . (2.13) 

It follows that v 1 = Av l — c v v 2 from which we deduce c v = v 2 ■ (Av 1 — v 1 ) = v 2 ■ (Q x v 1 ) — 
v 2 ■ v 1 = Q ■ v 3 + v 1 ■ v 2 . Thus we obtain the useful formula 

c v {9) = fi(0, oj9 + O ) • v 3 (9) + v\9) ■ v 2 (9). (2.14) 

In addition, since Tr[J 2 } = -2, it follows from (ETEty that c v = -(1/2) Tr[JV T {AV -V)}. 

Remarks: 

(1) Equation ()2.7j) is equivalent to a change of basis for spin whereby S is expressed as 
S = V(9)S = Siv 1 + S2V 2 + S3V 3 so that S is the spin in the rotating frame represented 

by the matrix V. Moreover S = c v (9)JS so that S{9) = exp (j J ' c v {9')d9'j 5(0). 

Thus S3 is constant and S precesses around (0, 0, 1) at a nonconstant rate c v {9). From 
(|2.14|) the rate c v is, as one would expect, just a combination of the projection of Q 
onto v 3 and the rate, v 1 • v 2 , of rotation of v 1 and v 2 around v 3 . In our discussion of 
spin tune in later sections it will be useful to define a frame in which the spin precesses 
uniformly. From ()2.11|) we have 

S(9) = V{9) exp J {c v {9') - v)dS^j exp (JvO) V T (0)S(0) , 

where v is an arbitrary constant. Thus with the change of basis 

S = U{9)S , U{9) = V{9) exp J (c v (9') - v)d9' 

we obtain S = vJS and we have defined a frame U whose third column is a spin 
trajectory and in which spin has a constant precession rate v. In the following we will 

T 

be interested in the case where the mean c v = limT_ KX) (l/T) f c v {9)d9 of c v exists 



and we will choose v = v v where v v = c v mod 1 and is in [0,1). Thus we can write 
c v W = (^) + v v + w h ere c v represents the fluctuating part of c v with zero mean 
and where the integer k v is chosen such that v v e [0, 1). 
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(2) The ideas in Remark 1 lead to some precise definitions. On a given torus, let U(9) G 
SO (3) be such that the principal solution matrix at 0o can be written as 

$(0; (f> ) = U{9) exp(Jv6)U T (0) , (2.15) 

where v is constant and in [0, 1). Then U is called a uniform precession frame (UPF) 
at <fio and u, which is uniquely determined by U, is called the uniform precession rate 
(UPR) for U and is denoted by u s (U). We then call 1)2.15)1 a standard form of the 
principal solution matrix. 

Under certain conditions which will be described in Definition 5.5, u s {U) will be called 
a spin tune. Note that, due to ()2.15j) . U(9) is smooth in 9 and satisfies the ordinary 
differential equation: 

U(e) = A{9-MU{0)-vU{0)J, (2.16) 

where v = v s (U). In particular by ()2.16)) the vector described by the third column 
of U obeys (|2.1|) so that it is a spin trajectory with unit length. Moreover, for every 
constant v the initial value problem defined by (J2.16)) and the arbitrary initial matrix 
17(0) G 50(3) has the unique solution U{9) = $(9; <f> o )U(0) exp(-Jv6). Thus every 
solution U{9) of (l2~T6l) with £7(0) G 50(3) and v G [0, 1) is a UPF at <p and its UPR 
y s {U) equals u, i.e. 

U(9) = A{9- <j) )U(9) - u s (U)U(9)J . (2.17) 

Because Tr[J 2 ] = —2, for every UPF U one has the useful formula 

u s (U) = —(1/2) Tr[J (U T AU — U T U)} , (2.18) 

which follows from (|2.17|) . Note that the interval [0, 1) is just a matter of choice - any 
convenient half open interval of length 1 could be chosen, e.g. (0, 1]. 

(3) In this section and in the rest of this paper, the concepts of orthonormal reference frame 
and 50(3) matrix are interchangeable. Moreover the elements of the columns of such 
a matrix are just the components of the unit coordinate vectors of the corresponding 
frame, as for example in ([2.6)1 . Thus we will often identify the columns with such 
vectors. 

(4) It can be shown that if V : R — > 50(3) with V =: [v 1 , v 2 , v 3 ], is a smooth and 27r- 
periodic function (e.g. a 27r-periodic UPF) then the fractional part of 

X '■= (1/2tt) Jg 27r d9v l {9) ■ v 2 (9), is independent of v 1 and t> 2 , i.e. the fractional part of 
X only depends on v 3 . If v 3 were represented in the "spinor formalism" ^3] then it 
would be found that the fractional part of x is the geometrical phase of v 3 in the sense 
of [Ml. 

Finally, we summarize the basic eigenstructure of an 50(3) matrix R and its exponential 
representation in terms of a skew-symmetric matrix B in the following lemma. 
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Lemma 2.1 a) Let R be a 3 x 3 matrix in 50(3). Then a real number \x G [0, 2tt) and a 
SO (3) matrix W exist such that the spectrum of R, X(R), is the set {e lfl , e~ v , 1} (whence 
the eigenvalues are on the unit circle) and such that 

(cos [A — sin /i \ 
sin/4 cos/i =Wexp(J», (2.19) 
1 J 

where in the second equality we used \2.1ty) . Furthermore, any such W = [w 1 ,w 2 ,w 3 ] 
satisfies the relations R^w 1 ± iw 2 ) = exp(=Fi//)(w 1 ± iw 2 ) and Rw 3 = w 3 . Also, by \2.19i) . 
R = e Bfl where the matrix B = WJW T is skew-symmetric. 

b) Conversely, if B is a real skew-symmetric matrix and if \{B) = {i, — i,0}, then a SO (3) 
matrix W exists such that B = WJW T . 

The proof is elementary. See for example [SZ1 dB] • Note that the relation Rw 3 = w 3 simply 
means that w 3 lies along the "axis of rotation" for R. 



3 Spin Motion on the Closed Orbit 

In this section we consider the case J = 0, so that the A of ()2.2j) has no </>o or u dependence 
and is 27r-periodic and smooth in 9. This case corresponds to the 27r-periodic motion of the 
particle on the closed orbit and the following theorem applies. 

Floquet Theorem: For J = 0, there exists a v G [0, 1) and W G SO (3) such that 
the principal solution matrix defined by (|2.3|) and (J2.4|) is independent of 0o an d can be 
decomposed as 

**k <S>(6) = p(6)W exp(Ji>6)W T , (3.1) 

where the matrix p(6) G 50(3) is 2-K-periodic and smooth and where p(0) = I. Moreover, 
A($(2tt)) = {e ii/27T , e~ ii>2w , 1}. 

Proof: Since, at J = 0, A(9, 0) is independent of 4>, $(#; 0o) is independent of 4>o- From 
Lemma 2.1a, we know that there exist \i G [0, 2tt) and W G 50(3) such that A($(27r)) = 
{e*M e -i/, 5 X } and $ ( 2 7t) = where B = WJW T . Then with fi = 2vrz>, v G [0,1), 
A($(27r)) = {e ii>2n , e~ ii>2n , 1} and $(27r) = e m2n . A key property of the principal solution 
matrix is that 

$(# + 2tt) = $(fl)$(27r) , (3.2) 

which we can see by noting that the l.h.s. of ()3.2j) is a solution matrix of ()2.3j) by the 2n- 
periodicity of A and the r.h.s. is a solution matrix of ()2.3j) since $ is. They are equal, by the 
uniqueness of solutions to the initial value problem for ()2.3|) and ()2.4j) . since they are equal 
at e = 0. 
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Define p by p(8) := <$>{9)e~ me . Clearly p(0) = / and p{9) G 50(3) since $ G 50(3) and 
B is skew-symmetric. The periodicity of p is clear because 

p{e + 2tt) = $(0 + 2tt) exp(-Bv(6 + 2tt)) = $(#)$(2tt) exp(-5z>27r) exp(-£i>0) = p(0) , 

where (13 .2|) is used at the second equality. □ 

Using the Floquet theorem we now make several remarks concerning the principal solution 
matrix defined by ()2.3|) . (j2.4|) when J = 0, i.e. when A is independent of <po and 27r-periodic 
in 9. The " * " symbol which was specific to the theorem is not needed in the following 
remarks. 

Remarks: 

(1) If = p(9) exp(Bud), where v G [0,1), where the real skew-symmetric matrix B 
has the spectrum X(B) = {i, —i, 0} and where the matrix p{9) is 27r-periodic, then 
p, B, v will be called Floquet parameters. In particular v is called a Floquet frequency. 
Thus the Floquet theorem states that Floquet parameters exist and it implies that 
the principal solution matrix depends on two frequencies where the Floquet frequency 
emerges in addition to the circulation tune lu c — 1. Note that the Floquet parameter 
p{9) is a smooth element of 5*0(3) with p(0) = I. 

(2) If p,B,u are Floquet parameters as defined in Remark 1, then from Lemma 2.1b a 
W G SO(3) exists such that WJW T = B, whence $(#) = p(9)We Jud (p(0)W) T . Thus 
at every <p , p{9)W is a 27r-periodic UPF with UPR z/, as defined in Section 2. We 
conclude that every Floquet frequency is a UPR of a 27r-periodic UPF and that (recall 
Remark 2 of Section 2) a 27r-periodic unit-length function of 9 exists, which is a spin 
trajectory at every <f)Q. This is the uq mentioned in the Introduction. 

Conversely, if U is a 27r-periodic UPF, then by ()2.15|) 

$(#) = U{9) exp(Ju s (U)9)U T (0) =: p(9) exp(Bu9) , (3.3) 

where p,B,v, defined by 

p(0) : = U(9)U T (0) , B := U(0)JU T (0) , v = v s (U) , 

fulfill all conditions of Floquet parameters. Thus the UPR of every 27r-periodic UPF is 
a Floquet frequency, i.e. for 27r-periodic UPFs the UPR emerges as a Floquet frequency 
and thus as an extra frequency of the system. We conclude that the set of Floquet 
frequencies is identical with the set of UPRs which correspond to 27r-periodic UPFs. 

(3) To study the set of Floquet frequencies in more detail, we first consider two sets p, B, v 
and p, B, v of Floquet parameters, i.e. 

$(0) = p (6) exp(Bu9) = p(9) exp(Bi)9) . (3.4) 

From (dOJ at 9 = 2tt, we obtain $(2vr) = exp(2irBu) = exp(2?r J Bi>), so that A($(2vr)) = 
{e 2mv , e~ 2mu , 1} = {e 2nlu , e~ 2mv , 1}. Thus the set of Floquet frequencies has at most 



14 



two elements and (due to the Floquet theorem) at least one element. In particular, 
the set of Floquet frequencies has either one element (which then is equal to 0) or it 
has two elements u,l — u, both of them positive. Note that a Floquet frequency which 
is in [0,1/2] always exists. Moreover, as we will see when we introduce the general 
concept of spin tune with Definition 5.5, the set of spin tunes contains this set of Flo- 
quet frequencies (see Remark 4 of Section 6). This in turn allows us to select one of 
the Floquet frequencies as a "preferred" spin tune which we denote by u . In fact this 
corresponds to the customary choice in which e tu ° 2n is an eigenvalue of <I>(27r) [T41 I15j . 

(4) By the definition of Floquet parameters, v = is a Floquet frequency iff $(27r) = /. 
From Remark 3 it follows that if is a Floquet frequency, then it is the only Floquet 
frequency. Furthermore all solutions of ()2.1)1 are 27r-periodic iff v = is a Floquet 
frequency. It follows by Remark 2 that for the case v = every 27r-periodic UPF has 
a zero UPR. See also Remark 4 in Section 6 and Theorem 6.4. 

The Floquet theorem does not give a method for constructing the principal solution 
matrix at J = since p is defined by $, W and v — it is a theorem about its properties. 
We now construct, by the method pointed out in Remark 1 of Section 2, a standard form 
for the principal solution matrix in order to compute Floquet parameters used in, or implied 
by, other treatments of spin motion on the closed orbit [To"! |2H EH] • We first show that the 
matrix (|2.fij) exists and can be chosen to be 27r-periodic in 9. 

Let £ 3 with length 1 denote the eigenvector for eigenvalue 1 of the 5*0(3) matrix <3>(27r). 
Then the solution of J2HJ) with 5 = £ 3 is 27r-periodic since S(9) = $(#)£ 3 = $(#)$(27r)£ 3 = 
$(6> + 27r)£ 3 = S(9 + 2ir). We now define v 3 (9) := $(0)£ 3 . To complete the construction of 
the matrix V we must now construct the 27r-periodic vectors v 1 and v 2 . To do this we can, 
for example, assume that a constant unit vector e exists such that v 3 (9) x e has no zeros. 
Then we can define v 1 (9) := v 3 (9) x e/\v 3 (9) x e| and v 2 (9) := v 3 (9) x v 1 (9). In any case 
the Floquet theorem ensures the existence of a smooth and 27r-periodic 50(3) matrix V(9), 
whose third column is v 3 (9): just choose a 50(3) matrix V(0) whose third column is £ 3 and 
then choose V{9) =p(9)V(0). 

We now consider an arbitrary smooth 27r-periodic 50(3) matrix V whose third column 
is a solution of ()2.1)1 and we construct a 27r-periodic UPF and the corresponding standard 
form for the principal solution matrix. From (|2.9jl and (|2.1()jl c v (9) is 27r-periodic and can 
be written as c v = c v + c v where, as before, c v and c v are the mean and zero-mean parts 
of c v . From Remark 1 of Section 2 we can write 

U{9) = V{9) exp(J(a(6) + k v 9)) (3.5) 

where da/d9 = c v and where v v e [0, 1) and the integer k v are chosen such that v v + k v = c v . 
The principal solution matrix is then given by ()3.3|) with v s (U) = v v . Because the mean of 
c v vanishes, the integral L c v {9')d9' is 27r-periodic in 9 and thus exp(j7"a(#)) and U{9) are 
27r-periodic in 9. 

The standard form for the principal solution matrix that we promised is given by the 
first equality in ()3.3|) . where U is defined by (|3.5j) . and the UPR corresponding to U is given 
by u s (U) = c v — k v . In ()3.5|) we have a procedure to calculate U{9) and u s (U) and thus 
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knowing a periodic solution of (|2.1|) . From £/"(#) and ^ S (Z7) one can construct Floquet 
parameters as in Remark 2. 

Remarks: 

(5) From Remark 2 it is clear that the unit-length 27r-periodic spin trajectory n Q exists 
and from Remark 4 it is clear that the direction of n is not unique iff vq = 0. 

(6) The above construction for v 3 (9) is easily done numerically. First integrate (|2.3J) from 
9 = to 9 = 2tt numerically to determine $(27r) and then solve the linear system 
$(27r)£ 3 = £ 3 . Then v 3 can be found on a grid of points in [0, 2n] by numerically 
integrating (j2.1|) with Sq = £ 3 . This is the way that n (9) is constructed in SLIM 
and other related spin codes. 

(7) Since v 3 is a 27r-periodic solution of ()2.1|) so is — v 3 . One can therefore replace V = 
[v , v 2 , v 3 ] by V = [v , — v 2 , — v 3 ]. By (12.14)1 it follows that with this replacement 
c v becomes — c v . Thus if [v 1 , v 2 , v 3 ] leads to u v G (|, 1) then [u 1 , — 1> 2 , — 1> 3 ] leads 
to v v G (0, =), so that one can choose V in (|3.5jl so as to put the UPR of U in [0, |]. 



4 Quasiperiodicity and a Diophantine Condition 

In the previous section we set up the Floquet form for the principal solution matrix for the 
case J = where A(9) is 27r-periodic. However, one of our aims is to obtain an analogous 
form when J ^ 0, i.e. when A contains the frequency vector u in addition to the circulation 
tune 1. For this and other purposes we need to introduce the concept of quasiperiodicity. 

A periodic function has one basic frequency whereas a quasiperiodic function has a finite 
number of basic frequencies (tunes) denoted by v = (z/ 1; u 2 , . . . , z/ fe ). We give the following 
definition. 

Definition 4.1: a) A function / : R — ► R is said to be quasiperiodic with tune vector v in 
M. k if a continuous and 27r-periodic function F : R fc — > WL exists and 

■k-kk f(9) = F(u9) . (4.1) 

A function / : R — > C is said to be quasiperiodic with tune vector v in M fc if its real and 
imaginary parts are quasiperiodic with tune vector v in R fc . A real or complex matrix valued 
function is said to be quasiperiodic with tune vector v in M. k if its components are. 

b) The spectrum of a quasiperiodic / is defined by A(/) := {A G R : a(f, A) ^ 0}, where 

a(f, A) := lim )- f f{9) exp{-iX9)d0 . (4.2) 

c) The mean of / is defined by 

/:=a(/,0), (4.3) 
and the zero-mean part of / is denoted by / and defined by / =: / + /. □ 

The class of functions so defined will be denoted by Q(u; k), where we refer to u as the 
tune vector and k as the order. If either v or k are obvious from the context or not relevant 
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we may omit either or both. We prove that a(f, A) exists in Lemma 4.3d. Note that for k = 1 
all functions in Q are periodic. A tune vector v is said to be nonresonant if the equation 

m ■ v :— m\V\ + • • • + m^v^ = 0, 

where m G Z fc (i.e. m is a k vector of integers), has m = as the only solution. If there are 
nontrivial solutions then v is said to be resonant. For some applications, e.g. as in Lemma 
4.3d, a tune vector may be assumed to be nonresonant, since otherwise the order can be 
reduced, as the following lemma shows. 



Lemma 4.2 If f is in Q(y\ k), then a nonresonant v exists such that f is in Q(p\ k), where 
k<k. 

Proof: Let / G Q(v; k) and let / : R -»• R. Then f(B) = F(u6) where F : R k -> R 
is continuous and 27r-periodic. If v is nonresonant, then nothing has to be proved. If 
v = 0, then / is constant, so that / G Q(l; 1). If v is resonant and nonzero then k > 2 
and there exists nonzero m G Z fc such that m ■ u = 0. If 7^ 0, then F : R fc_1 — > R, 
defined by ^-i) := F(— m^Zi, — m^Zk-i, ra\Z\ + ... + rrik-iZk-i), is continuous 

and 27r-periodic and f{6) = F{v9) = F{i>9) where v := — (l/m&)(i/i, ^-i). The case 
= can be dealt with analogously because m is nonzero. We have therefore shown that 
/ G Q{p] k — 1), where v 7^ 0. If z> is resonant, we can repeat the above procedure until 
a nonresonant tune vector is obtained. For complex valued functions and matrix valued 
functions the proof proceeds analogously. □ 

A different proof of Lemma 4.2 can be found in j^OJ Appendix 3.8]. 

The order in every Q can also be increased since the F in ()4.1j) can be viewed as a function 
with domain Thus for an arbitrary real vector \x in R , Q(v; k) C Q(u, \i\ k + I). For 

later use we note that if a is an integer, exp(J~a8) is 27r-periodic in 6 and therefore in Q(l; 1). 

The class of functions Q(v) obviously depends on u, but two different v can generate the 
same class. For example, Q is unaltered if v is replaced by v = Qv where Q is a matrix 
with integer entries and determinant ±1, i.e. / G Q(u;k) iff / G Q(v,k). To see this let 
/ G Q(v; k). Then there exists a function F such that (|4.1j) holds. Let H(z) := -F(Q _1 2;). 
Then H is 27r-periodic since Q~ l has integer entries, and / G Q(0; k) since f(9) = 11(06). 
The converse is proved in the same way. 

The periodic function F in Definition 4.1 has the Fourier series 



E 



F m exp(im ■ z), (4.4) 



where the Fourier coefficients of F are defined by 

I /*2-7r rift 

F m := 77— rr / • • • / F(z) exp(-im • z)dz x ■■■dz k (4.5) 
(27r) fe Jo J 

and where m • z := mi^i + • • • + m^Zk. If F is of class C k and if we define 

Sn(z) ■■= ^2 F rn,exp(im ■ z) , (4.6) 

mez fc 

||m||<JV 
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then the sequence {SW} converges uniformly to F on M fe [JTJ p. 411]. Furthermore, the 

sequence {/jv}, defined by /jv(#) := S^{u9), converges uniformly to /. Thus we have the 
representation 

f{9) = lim f N (9) = lim V F m exp(z(m-z/)£), (4.7) 

A/ — >rvi /V — >r*n ' * 



N^oo N—fOO 

m€Z. k 

||m||<JV 



for / in terms of the Fourier coefficients of F. Here the norm ||m|| of the integer vector m 
is the max norm, i.e. | \m\ | := max{\m,i |, |rafc|}. 

For our definition of quasiperiodicity F is only required to be in C°, not in C k . Thus in 
general (|4.7J1 does not apply and we need the following lemma. 



Lemma 4.3 a) Let F be a continuous and Ik -periodic real or complex matrix valued func- 
tion on R fe . Then the sequence {£tv} ; defined by 

E N (z) := ^2 A N,m F m exp(im- z) , (4.8) 

7 

l< 

converges uniformly to F on~R k , where 



||m||<iV 



yj iV + 1— \m n \ 

AN > m '- 11 w+i 

n=l 



b) If f G Q(i / ; fc) and if F is a continuous and 2ir-periodic function on M. k such that f(6) = 
F{v9), then 

f{6)= lim EjvO/0) = lim V A Nm F m exp(i(m -u)9), (4.9) 

N^oo N—^oo ' * 

mSZ fe 
|m||<JV 

where the convergence is uniform. Thus ()4.9|) generalizes ()4.7|) to this less smooth case. 

c) Let v be nonresonant under the conditions of Lemma 4-3b. Then the mean of f , defined 
by ([4.3)) . exists and is given by 



1 



i-2-K r 2-rr 



f = F ° : =J^AkJ o •••] F(z)d Zl ---dz k (4.10) 
and the Fourier coefficients F m of F in ()4.5j) satisfy the relation 

F m = lim ^ / f{9)e- i{m v)e d9= lim }- [ F{v9)e- i{m v)e d9 . (4.11) 

In particular, f = implies F = 0. 

dp Lei / 6e m fc). TTien a(f, A) exzsfo /or a// A. 

e) Let f be in Q{v] k). Then A(/) C {m ■ v : m G Z fc }. 
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Proof of Lemma 4- 3a: See for example [42]. □ 
Proof of Lemma 4-3b: A simple consequence of Lemma 4.3a. □ 
Proof of Lemma 4 -3c: That / exists and equals Fq is proved for example in O Section 
10.3]. The proof of 1)4.11)1 is similar. □ 
Proof of Lemma 4-3d: Let / G Q.{y\ k). Then by Lemma 4.2, /exp(— iX ■) G Q(v] k), where v 
is nonresonant and where / exp(— iX ■) denotes the function which maps 9 to f{9) exp(— iX 9). 
Applying Lemma 4.3c we conclude that a(f, X) exists for all A. □ 
Proof of Lemma 4-3e: Since a(Ejv(z/ •), A) and a(f, X) exist by Lemma 4.3d we have 

\a(E N (u •), A) - a(f, A)| = \a(E N (v) -f,X)\ 

= \ lim )- [ (Z N (v9)-f(9))exp(-tX9)d9\ < sup \Z N (v9) - f(9)\ , (4.12) 

where we note that £tv(z/ ■) and / are bounded functions and where | • | denotes the Euclidean 
norm. It follows that 

lim a(Z N (u-),X) = a(f,X) , (4.13) 

since the convergence in (J4.9j) is uniform. If A G" {m ■ v : m G Z k }, then a(Ejv(z/-), A) = 
and g3S|) gives a(/, A) = 0. Thus if A G" {m ■ v : m G Z fe }, then A G" A(/), so that 
A(/) C {m ■ v : m G Z k }. □ 

Remarks: 

(1) Lemma 4.3a is a multidimensional generalization of Fejer's theorem and it shows that 
a continuous function can be recovered from its Fourier coefficients. For k — 1, £/v m 
(14. 8 j) is the Cesaro sum of the sequence {Sat}. For k > 1 it is a natural generalization 
of Cesaro summability and is one of many techniques for summation. The claim in 
Lemma 4.3c that / = Fq, is a "flow" version of Weyl's equidistribution theorem which 
is one of the main theorems of ergodic theory 26J, 27, Chapter 3],[43j. Of course, if F 
is in C k then Sn and E^r converge to the same function F. 

(2) Definition 4.1 gives / in terms of F, and Lemma 4.3a and (j4.11|) show how to recover 
F given / with nonresonant v. Furthermore Lemma 4.3b gives a representation for / 
even when ()4.7|) is not valid. If v is resonant then Lemma 4.3 does not allow F to be 
recovered. However we can find F for the nonresonant tune vector v of Lemma 4.2. 

(3) If / is quasiperiodic, then by Theorem 4.3d / exists, whence /:=/ — / is quasiperiodic 
and the mean of / exists. In particular / = / — / = which justifies calling / the 
zero-mean part of /. 

(4) From (|4.8jl and (|4.11j) we see that F, for nonresonant u, is determined by its values on 
the curve D' := \v9 : 9 G R}. That this is to be expected can be shown as follows. 
Because F is 27r-periodic, the knowledge of F at the points of D' implies that F is 
even known at the points of the larger set D := {u9 + 2nM : 9 G E, M G Z k }. But D 
is dense in M fc since v is nonresonant and any continuous function is determined by its 
values on a dense set. 
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The following definition will be convenient in some of our proofs: 



Definition 4.4: A function / : R — ► R is said to be almost periodic if a complex sequence /„ 
and a real sequence A n exist such that the sequence of functions X^=i fn exp(iA n 6*) converges 
uniformly to / as N — > oo. A function / : R — > C is said to be almost periodic if its real and 
imaginary parts are. A real or complex matrix valued function is said to be almost periodic 
if it is almost periodic in each component. □ 

Remark: 

(5) Due to Lemma 4.3b, the quasiperiodic functions form a subset of the almost periodic 
functions. In contrast to the quasiperiodic functions, the almost periodic functions 
generally have no finite set of basic frequencies. The Fourier representation of such a 

oo 

function is given by f n exp(i\ n 9), where {A n } is the set of frequencies. A quasiperi- 

71=1 

odic function can be written in this way, but ()4.7|) is often more convenient for analysis. 
Almost periodic functions were introduced by H. Bohr and a standard modern treat- 
ment can be found in jj^, and in summarized form in the appendix of [3*5] . 

We have now provided the basic machinery needed to deal with the quasiperiodic func- 
tions appearing in the following sections. In the remainder of this section we treat a special 
aspect of quasiperiodicity, namely the so called small divisor problem and we show how one 
solves this problem by using sufficiently differentiable functions. This material will only be 
needed in Theorems 6.5c-d. 

Integration as in the definition of a(9) in Section 3 (see (|3.5Jl ) is, of course, a smoothing 
operation and the integral of a periodic function with zero mean is periodic. It is perhaps 
surprising then that the integral of a quasiperiodic function with zero mean is not necessarily 
quasiperiodic. This is an example of the so-called small divisor problem. To see this let v 
be nonresonant and let / be in Qiy\ k) with zero mean with its associated F in C k . Define 
g{9) = f Q f(8') d& for / given by ([4.7)1 . Because of the uniform convergence in (|4.7|) we can 
interchange limit and sum giving 

g(9) = lim V -i(m ■ z/) _1 F m ( exp(i(m • v)6) - 1 ] . (4.14) 

0<||m||<JV 

However, although g is well defined, it may not be quasiperiodic if k > 2. The source of the 
problem is that the divisor m ■ u, while not zero, can be arbitrarily small and this can lead 
to an unbounded g, thereby contradicting (4.1). We give such an example in Section 8. To 
find a sufficient condition for quasiperiodicity we define 

G N (z) := ~^( m ' u )~ lF ml exp(im ■ z) - 1 J , (4.15) 



mgZ fe 
0<||m||<JV 



and note that g{9) = UmN^ooGN^O) so that Gn converges on the lines z = uO + 27rZ in 
M fc where I is an integer. Now suppose that G^ converges pointwise in R fc to a function G. 
Then G is 27r-periodic and g{6) = G(u6). Since g is continuous, G is continuous along the 
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lines z = vQ + 2nl, but Definition 4.1 requires that G be continuous everywhere on R fc . So 
we need a condition that ensures both pointwise convergence of Gjy and the continuity of 
the limit. Uniform convergence guarantees both, and a sufficient condition for this is a so 
called "Diophantine condition" which places a condition on how close m ■ v can be to zero. 

In the following we will be concerned with the tune vector v = (1, o>) where 1 is the 
circulation tune and uj G M. d is the orbital tune. Thus we formulate the Diophantine condition 
and the convergence of in terms of v — (1, to) and k — d+ 1. See [JUJ Appendix 4] for 
some details on Diophantine conditions and [13] for further discussion and its use in another 
beam dynamics context. 

Definition 4.5: The tune vector uj G M. d is said to satisfy a Diophantine condition if, for 
positive r, uj is in the Diophantine set fl(r) := U 7 e(o l] ^( r '7) "where f2(r, 7) is the set 

{w G M d : |m • > 7||m||- T ,m G Z d+1 , m + 0} 

□ 

Note that the symbol f2 is also used in the Introduction for the precession vector. However 
its meaning should be clear from the context. 

It follows from the definition that fl c (r, 7) = Umez^+^o} ^( r > 7> m )> wnere ^( r > 7> m ) : = 
{w G l 11 : |m ■ (l,w)| < 7||m||~ T } • For fixed (r, 7) and for each m the "resonant zone" 
Z(t, 7, m) is either empty or a thickened d — 1 dimensional plane centered on the resonant 
plane, m ■ (1, o>) = tuq + miCo>i + ■ • ■ + vn^iOd = 0, with thickness proportional to 7. For 
example, when d = 1, ||m|| = 1 and 7^ the corresponding zones are intervals centered 
on the three points in 1R: lo\ — r\ where 77 G { — 1,0,1}, each with thickness 27. When 
d = 2, \\m\\ = 1 and (mi,m2) 7^ the corresponding zones are centered on the twelve lines 
in M 2 : uji = 77, CU2 = r], cj 2 = ±^1 + f] where r] G { — 1, 0, 1} with thickness either 27 or -v/27. 
More generally, if m = (mo, m) with tuq G Z and rh G (Z d \ {0}) then one can show by using 
rotations in M. d that Z(r, 7,m) can be rotated into the set {uj G K d : |g + u;i| < ^ T }. 
Thus the thickened d — 1 dimensional plane 2(r, 7,m) has thickness 27 1 |m| | _r /|m|. For m 
such that m = the resonance condition cannot be satisfied and we have Z(t, 7, m) = 0. 
Note also that Z(t, 7, m) is undefined for m = 0. 

Definition 4.6 (orbital resonance): We say that the torus at J is off orbital resonance if 
(l,u(J)) is nonresonant (otherwise we say that it is on orbital resonance) and this is certainly 
the case if oj(J) G Vt{r). Usually, (jl.5j) is said to be resonant if uj is resonant. Our usage is 
different because our basic system is (|2.1j) . which includes the circulation tune, 1. □ 

We can now interpret f2(r, 7) as the closed set in R d constructed by successively remov- 
ing the open resonance zones, corresponding to the resonance planes m ■ {1,uj) = 0, with 
increasing ||m||. Thus its construction is similar to the construction of a Cantor set. The 
resonance planes are dense in M. d and thus f2(r, 7) is small in the sense that it has an empty 
interior. However, we will show in the proof of Lemma 4.8 that it is large in the sense that 
for t > d + 1 the Lebesgue measure of its complement relative to R d is proportional to 7 (in 
the sense of (J4.23|0 which can be arbitrarily small. We could take our Diophantine set to be 
Q(t, 7) for 7 small as in [IS]. Here we take the larger set Q(t). 
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Now ii u E fi(r, 7) then |(m • (l,u;)) _1 F m \ < 7 _1 | \m\ \ T \F m \ and thus Gat converges 
uniformly if ||m|| T |i^ m | converges. Thus the Diophantine condition leads to a simple 
sufficient condition for the quasiperiodicity of g. In this context we now state and prove the 
following lemma which addresses the differentiability of F . 



Lemma 4.7 Let F : R d+1 — > R be of class C n and 2ir -periodic and let F = 0. Let to E Q(t) 
where < r < n — d — 2. Then Gn '■ M. d+1 — > IR, given by (j4.15J) with v = (1,&>), converges 
uniformly on M d+1 to a smooth function G which is 2tt -periodic. Moreover, VG(z) • = 
F{z). 

Proof: Because a constant M > exists such that ||m|| n |F m | < M (see |4*T] p. 409]), we have, 
for every r > 0, 

TV 

\\m\\ T \F m \<M ^ \\ m \\ T - n = M^2f- n ^ 1 - ( 4 - 16 ) 

:\\m\\<N 

A combinatorial argument gives 



m£Z d + 1 m6Z ti + 1 j=l 

0<\\m\\<N 0<||m||<JV 



1 <3 d+1 (d+l)f . (4.17) 



l»»ll=J 



Moreover, if uj E f2(r, 7), then |(m ■ (1, *| < 7 1 ||m||' r and ()4.16|) and (|4.17jl give 



/v 



£ |(m-(l,c))~ 1 F m |<3 d+1 (rf+l)M 7 - 1 ^r-" +d . (4.18) 



0<||m||<iV 



Since t < n — d — 2 we conclude that the l.h.s. of 1)4.18)1 converges as N — > 00 for every 
w G O(r). It follows that Gn converges uniformly to a continuous function G which is 
27r-per iodic. 
From (J4.15)) we have 

S/G N (z) = Y ( m ' (l,^))" 1 i r mexp(im • z) m . (4.19) 

m£Z<2+l 
0<||m||<JV 

Then if uj E f2(r, 7) and repeating the above argumentation we find 

\VG N (z)\< Yl (l,u))- l F m \ \m\ < Vd+1 l( m- 0->v))~ lF m\ INI 

mgzd+1 m6Z d + 1 
0<||m||<AT 0<||m||<JV 

TV 

< 3 d+1 (d + 1) 3 / 2 M 7 - 1 j T+1 ~ n+d ■ (4.20) 



22 



Thus, if uj G fi(r), VGV converges uniformly and a standard result (see jlHl 8.6.3]. |4T| 
p. 117]) means that G is C 1 and 

VG(*)=Iim > (m- {l,u)Y l F m exp(im • z) m . (4-21) 

N—*oo ' ^ 

0<\\m\\<N 

It follows that VG(z) • = limjv^oo VGjv(^) • (1,^0 = F(z) since, with n > d + 1, the 

Sn in f)4.6|) converge. □ 
Remark: 

(6) To prove ()4.17|) we define the sets s := {in G Z d+1 : ||m|| = j} and Sj := {m G s : 
\rrii\ = j}. It follows that s = [J^ 1 Sj and that Sj contains 2(2j + l) d elements. Then 
s contains no more than 2(d + 1) (2j + l) d elements and because 

2(2 j + l) d < 3(2j + l) d < 3(3j) d = 3 d+1 j d , 

we conclude that s contains no more than 3 d+1 (d + l)j d elements, thus proving ()4.17)) . 

Lemma 4.7 provides the basic framework that we need for discussing the uniform con- 
vergence of the sequence Gn- In particular it shows that as n increases beyond r + d + 2 
the small divisor problem loses much of its potency. This comes as no surprise because the 
inequality ||m|| n |F m | < M implies that the Fourier coefficients decrease with increasing 
||m|| more rapidly as n increases. Then with growing n the small divisor in ()4.15|) can come 
closer to zero without destroying the convergence. 

However, although Lemma 4.7 takes the mystery out of the working of the Diophantine 
condition, it puts the burden on determining which uj are in the set fi(r) and it is not so 
easy to decide, off orbital resonance, if uj is in Q(t). But some relief comes from the following 
lemma, which shows that if r > d+ 1 then the complement, fi c (r), of ^(t) is a small set in 
terms of Lebesgue measure. If in addition r < n — d—2, the sequence GV converges uniformly 
for almost every uj. For these two conditions to be consistent we thus need n > 2d + 3. 
These results will be central to the statement and proof of Theorems 6.5c-d. 

Lemma 4.8 If r > d + 1, then /z(f2 c (r)) = ; where \i denotes the Lebesgue measure. 

Proof: Let B(R) := {uj G R d : \uj\ < R} and define &(R,m) := Z(r,-f,m) H B(R). Then 
from Definition 4.5, fi c (r, 7) fl B(R) = Umez^Uo} ^(R, m )- Note that Z(r,j,m) (hence 
&(R,m)) is undefined for m = 0. 
We will show below that 

fx(&(R,m)) < 2~fR d - 1 a(d-l)\\m\\- T , (4.22) 

where a{n) := n n ^ 2 /T(n/2 + 1) (note that a(d) is the volume of B(l)). 
Assuming the validity of ()4.22j) . 



H m c (r, 7) n B(RU < 2 1 R d - 1 a(d - 1) ^ 



m£Z d + 1 \{0} 
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As before (see (Q7I) ) 



E inr = Er T E i<3 d+ v+i)E/~ r - 

meZ d+1 \{0} i=l mezd+i j=l 

l|m||=j 

Therefore 

io c (r )7 )nB(E)] < 2jR d ~ 1 a(d - l)3 d+1 (d + l)^/- T , (4.23) 

and for r > d + 1 the series converges. Now f2 c (T, 1/p) decreases monotonically to ^ c (t) as 
the positive integer p — > oo. Therefore 

/i^ c (r) n £(fl)J = lim /i^ c (r, n S(i2)J = 

by continuity and the finiteness of Lebesgue measure restricted to B(R) where the second 
equality follows from 1)4.23)1 . Since this is true for all R we obtain the required result. 

To prove (j4.22j) we first recall that Z(r, 7, m) is empty if m = (m , 0, 0). Thus &(R, m) 
can only be nonempty if m = (mo, rh) with mo G Z and m G (Z d \ {0}) so that we only have 
to consider this case. Using the fact that volumes are invariant under rotations it follows 
from the remarks after Definition 4.5 that 

u(&(R,m)) <^({uj G B(R) : + uj x \ < 7||mir r } ) , (4.24) 

V \ m \ J 

where we also used the fact that l/|m| < 1. If ci = 1 then ()4.24|) and the fact that a(0) = 1 
immediately yield ([4.22)1 . Furthermore, for d > 2, we conclude from (|4.24j) that 

li(e(R,m)) < fi({cu G R d : yju^ + ... + < fl, < 7||m||- T }^ 

d . I , ,2 1 1 , ,2 D I, , I ^ „,| l™l l-r-i \ _ r> ud-l 



= fMl{u G R d : y/u$ + ... +u 2 d < R, {u^ < 7||m||~ r }J = 27#*" A a(d- l)||m||" T , 

where the r.h.s. is just the volume of the cylinder in M. d with height 27||m||~ r and radius R 
and where in the first equation we used the fact that volumes are invariant under translations. 
This completes the proof. □ 

5 Introduction to the Quasiperiodicity of the Spin Mo- 
tion. 

Definition of the Spin Tune 

We can now continue the study of the principal solution matrix <&(9; 4>o,oS) for (j2.1j) which 
is defined by the initial value problem 

— = A(0;0 o ,u;)$, $(O;0o,w) = J (5-1) 
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where 



A(9;4>o,uj) :=A(e,cu9 + 4> ) . (5.2) 

In the language of Section 4 A(-; 4> , ui) is quasiperiodic in Q(l, u>; d + 1). One of the aims of 
this paper is to define the spin tune for this system. We are guided by the special case J = 
of Section 3 and therefore, and as will become clear below, our emphasis in this section is on 
principal solution matrices, where $(-;0o,^) is quasiperiodic and can be written in the 
standard form, 

**k $(0; 0o, a) = U(6) exp(Jv6)U T (0) , (5.3) 

where U is a quasiperiodic 5*0(3) matrix in Q(l, ui;d + 1) and v G [0, 1) and both may 
depend on (<pQ,u). In that case all solutions of (j2.1J) are in Q(l, u), v\ d + 2) with a very 
simple frequency structure in the tune v. Such principal solution matrices fulfill a generalized 
Floquet theorem. 

There has been extensive study of the equation 8^/89 = A(6;y)§, where the matrix 
A(-;y) is almost periodic and y is a vector of parameters, one goal being to find conditions 
under which almost periodic solutions exist (see for example [44, 48 ). However our problem 
()5.1|) is quite special because of the parameter dependence of A induced by A(9, 0) in ()5.2|) . 
In fact, for every integer iV we have 

A(6 + 2vriV; O , uj) = A{6; 2ttNu + O , w) , (5.4) 

which follows from ()5.2|) and the 27r-periodicity of A(8, 0) in 9. We will return to this in 
Section 9 where we will see that condition (|5.4|) has useful consequences for the spectrum 
A($(-; (f) ,uj)) of the principal solution matrix. 

The case J = where A(-;(p ,u) is 27r-periodic, i.e. in Q(l; 1) and independent of O 
and u, was discussed in Section 3. There we found a solution of (j2.1J) in Q(l; 1), namely 
v 3 (8) = p(0)C, 3 - All other linearly independent solutions are in <2(1, v\ 2). The latter follows 
directly from the Floquet theorem, but more importantly also from the construction of 
the 27r-periodic UPF using v 3 which led, with U in Q(l; 1), to the representation ()5.3|) . 
It is this construction that points the way for a generalization of the Floquet theorem to 
A e Q(l,w,d+1). 

We begin with the following proposition. 

Proposition 5.1 Consider the initial value problem 

S = A{6)S, 5(0) = So, (5.5) 

with a real skew-symmetric 3x3 matrix A 6 Q(v; k) and make the following assumptions: 

a) Equation (J5.5j) has a nonzero solution v 3 in Q{v\ k). 

b) There exists a smooth 50(3) matrix V = [V, t> 2 , t> 3 ] in Q(v; k). 

c) Let c := v 2 ■ (Av 1 — v 1 ). Since c G 0,{v\ k) its mean, c, and zero-mean part, c, exist (recall 

e 

Remark 3 of Section 4)- We assume that J c(8')dd' G Q(u;k). 

o 
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Then all solutions of ()5.5|) are in Q(v, Vk+x] k + 1) where Vk+i = c mod 1 and Vk+i £ [0, 1). 

Proof: From (c) there is an integer I such that c(6) = c + c(9) = Vk+i + I + c(9) and by the 
same calculation that leads to (j2.11j) the principal solution matrix for (|5.5J) is 

$(0) = V{9) exp J c(9')d9' + W)j exp {Ju k+l 9) V T (0) . 

Clearly, $(6>) G Q(i/, fjt+ii + 1) and thus all solutions are quasiperiodic as stated. □ 

However as pointed out in Section 4, even if conditions a) and b) are fulfilled condition 
c) might still be problematic since the integral of a zero-mean quasiperiodic function may 
not be quasiperiodic. 

To continue our discussion of the quasiperiodic structure of the principal solution matrix 
defined in (15. lj) and its potential representation in (J5.3)) we elaborate on the definition of 
UPF and UPR from Remark 2 in Section 2. 

Definition 5.2 (UPF and UPR): Consider ()5.1|) for a fixed O on a fixed torus J . If v G [0, 1) 
and U G 5*0(3) exist such that the principal solution matrix at 0o can be written as in 
(15. 3 j) . then the orthonormal coordinate system represented by U is called a UPF (uniform 
precession frame) at 0o arid v =: v s (U) is called the UPR (uniform precession rate) for U. 
If U G Q(l,u;d + 1), then it is called proper. A UPR is called proper if it corresponds 
to a proper UPF. The set of all proper UPFs at O is denoted by il(0o) and the set of all 
proper UPFs on the torus is given by il := U^eR"* ^(M- m addition we denote the set of 
all proper UPRs at 0o by H(0 O ) and then the set of all proper UPRs on the torus is given 

b y s : = LU )eK dH(0 o ). □ 

For further background material see Remarks 1 and 2 in Section 2. Recall that in Definition 
5.2 "on a fixed torus Jo" means that the orbital tune has the value u(Jq) and that the spin 
motion is characterized by the function A(-,-,Jo). Recall also that we often suppress the 
symbol Jo- From the discussion in Section 2 we know that UPFs always exist. However the 
existence of a proper UPF imposes additional constraints. For example Proposition 5.1 gives 
rise to the existence of a proper UPF for (J5.1j) (and a proper UPR) but the conditions are 
rather severe. If a proper UPF U at O exists, then clearly $(-;0 o ) £ Q(l> ^, ^(U); d + 2). 
However it is unknown if the converse holds although this is a reasonable conjecture. Note 
that for each proper UPR there is an infinite number of proper UPFs: if U is a proper UPF 
at 0o then U exp(J'y) is a proper UPF at (fio if y is a real constant. 

It is easy to see that if a proper UPF U exists at a given <po, then for all v G [0, 1) such 
that v = ±u s (U) +j +j ■ to, where (j , j) G Z x Z d , U(9) := [u 1 (9), ±u 2 (9), ±u 3 (9)} 
■ exp(— j7"(jo + j ■ &)9) is a proper UPF at that 4>o with u a (U) = v. This motivates the 
definition of an equivalence relation where v\ and v-i in [0, 1) are said to be equivalent - and 
we write V\ ~ z/ 2 - iff there exist (e,jo,j) G {—1, 1} x Zx Z d such that z/ 2 = £Vi+jo+j-uj. The 
equivalence relation partitions the interval into equivalence classes [u] := {// G [0, 1) : fi ~ u} 
such that [ui] = [1/3] iff v \ ~ v 2 and [ui] fl [z/ 2 ] = iff ^1 7^ z/ 2 . We note that if w has one 
irrational component then each [is] is a dense subset of [0, 1). To see this suppose that uj\ is 
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irrational. Then ev + j\OJ\ mod 1 is dense as j± varies over Z so that fi = ev + jo + j ■ oj is 
dense in [0, 1) as j an d j vary. 

From the above motivation for the definition of equivalence it is clear that if v\ G S(0o) 
then the equivalence v\ ~ v 2 implies v 2 G 2(</> ), i.e. [vi] C S(0 O ). Now suppose that U\ and 
U 2 are in il(0 o ) and that v s {U x ) = v x . Then exp{Jv s {U 2 )6) = U%{0)Ui{6)exp(Jv a {U{)0) 
■U?(0)U 2 (Q) so that thel.h.s. is in Q(l,to,v s (U 1 );d + 2). Thus i/ a (E7 2 ) = jo+j ■^+Jd+2^s(U 1 ) 
by Lemma 4.3e. It is plausible that — =tl- If that is so, then z/ x = v s (Ux) ~ ^(^2)- 
Then since C/2 is an arbitrary member of il(0o), E(^o) C [^1]. In fact this is the case, and 
the joint conditions [ui] C S(0 O ) and S(0 O ) C [ui] can be embodied in a theorem. 

Theorem 5.3 If S(0o) nonempty with an element v, then [v] = 5(0o)- 

Although all proper UPRs are equivalent at a given (fi on the torus Jo, it is not true in 
general that all proper UPRs on the torus are equivalent. But if they are, we then say that 
the torus is well-tuned and we call any v e 5 a spin tune. Since this situation is central to 
this paper, we delay the proof of Theorem 5.3 in order to formalize this definition. 

Definition 5.4 (well-tuned): A torus is said to be well-tuned iff the following two conditions 
hold: 

a) $(•; O ) has a proper UPF for each O £ i- e - each it(0o) is nonempty. 

b) Let U and U be in It, then u s (U) ~ v s (U). Thus for every v G 5, [u] D H and, by 
Theorem 5.3, [1/] C S so that then, [u] = S. 

A torus that is not well-tuned is called ill-tuned. □ 

Note that by Theorem 5.3 a torus is well-tuned iff the H(0o) have an element in common. 
This criterion is very convenient and we will use it for example in the proof of Theorem 6.3a. 
Note also that 5(0o) always contains at most countably many elements. In particular, if a 
torus is well-tuned then H contains at most countably many elements. 

Definition 5.5 (spin tune): Let a torus Jq be well-tuned, then each element of S is called 
a spin tune. Thus for each spin tune u, [u] = S. □ 

To prove Theorem 5.3 we need the following lemma. 

Lemma 5.6 Let exp(zco^), where Cq is a real constant, be in Q{v\ k) with k > 1. Then there 
exists n G Z k such that c = n ■ v. Moreover n is unique if v is nonresonant. 

Proof: First of all we observe that 

A(exp(ic ■)) = {co} • (5.6) 
Also, by Lemma 4.3e we see that A(exp(zco ■)) C {m ■ v : m G Z h }. Hence by fl5~6jl : 

c G {m ■ v : m G Z h } . (5.7) 
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Thus neZ' exists such that cq = n ■ u. Clearly n is unique if v is nonresonant. □ 
We can now prove Theorem 5.3. 

Proof of Theorem 5.3: Consider U =: [u 1 , u 2 , u 3 ] and U =: [u 1 , u 2 , u 3 ], both in iX((f>o), and 
define the smooth functions g± : R — > C by 



It follows from pTfjl that 



g±:= {u 1 ±iu 2 ) ■ {u 1 + iu 2 ) . (5.1 



g± = i[v s (U)±v s (U))g ± 



so that 

g±(6) = exp(i(v s (U) ± W(0) . (5.9) 

Because U and U are in Q(l, a;; d + 1), ()5.8|) ensures that g± are also in Q(l, u>;d+l). Now, 
if gf + is not the zero function then by (|5.9|) exp(i(v s (U) + v s (U))9) is in Q(l, cj; d + 1) and by 
Lemma 5.6 v s (U) + v s (U) = m ■ (1, cj) for some m G Thus either g + = or m G 

exist such that v s (U) + v s (U) = m- (1, to). Similarly either g_ = or m G Z d+1 exist such that 
v s {U) — u s (U) = m ■ By definition either g + or is different from the zero function, 

since otherwise u J ■ u k = for (j, = 1,2) which is obviously false. Thus m G Z a!+1 exist 
such that either u s (U) + v s (U) — m ■ (1, cu) or i/ s (C7) — ^ s (i7) = m • (1, uj). Hence e G {—1, 1} 
exists such that u s (U) = eu s {U) + m- i.e. S(0 O ) C [i/ a (?/)]. But by the remarks before 

Theorem 5.3 we also have that H(0 O ) 3 [^(^ )]■ n 

Remarks: 

(1) For the case of J = and arbitrary O eac h Floquet frequency u is in S(0 O ) (see 
Remark 2 in Section 3). We will see in Remark 4 of Section 6 that the torus J = is 
well-tuned so that every Floquet frequency vq is a spin tune. 

(2) As we will see in Remarks 13 and 14 in Section 6, matrices A(9, 0) can be found with 
which there are tori where S(0 O ) 7^ for every O but which are not well-tuned. 

(3) Consider a well-tuned torus. We say that the torus is on a "spin-orbit resonance" if 
G S. Thus on spin-orbit resonance the set of spin tunes reads as 

H = {u s G [0, 1) : uj s — jo + j • jo £ 2) j £ F° r ^ = 3 the spin-orbit resonance 

condition amounts to (|1.3|) . In general the condition takes the form 

co> s = m + mi w\ H +m d uj d (5.10) 

and the order of the resonance is |mi| + • ■ ■ + \m^. Thus on a well-tuned torus, if 
one 0o has resonant spin motion then all <f>o have resonant spin motion. Of course the 
same is true for nonresonant spin motion. This is a key aspect of being well-tuned: 
spin-orbit resonance is not defined in terms of a spin frequency u s that varies with 0o 
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in the sense that [u s ] is dependent on 0o- An important consequence of this definition 
of resonance is presented in Theorem 6.4. Clearly, the order of a resonance depends 
on the chosen u s and that in turn depends on the chosen UPF. Thus the order of a 
spin-orbit resonance is frame dependent. However, in practice the order is fixed by 
choosing a "preferred" spin tune and the corresponding UPF for every J. The way 
to do this will be explained in Section 10 where some numerical results on spin-orbit 
resonance are also mentioned. A numerical method for computing H is outlined in 
Section 9. 

If one component of uj is irrational and if the torus is well-tuned, then from the previous 
discussion the sets 5 and [0] are countable dense subsets of [0, 1). Thus even though the 
complement of S is uncountable and a number chosen at random from [0, 1) is unlikely 
to be a spin tune, S contains spin tunes arbitrarily close to any number in this interval. 
This simply implies that every spin tune (in particular the preferred spin tune) is close 
to a number in the resonance set [0] so that the system is always arbitrarily near to 
a high order resonance. Although the set [0] is a dense subset of [0, 1), the subset of 
low order spin-orbit resonances is not dense in [0, 1). Thus the closeness of a preferred 
spin tune u s to a number in the set [0] does not mean that u s is near to a low order 
spin-orbit resonance. 

The set Sj of proper UPRs on the torus J is a function of J which is uniquely deter- 
mined by the J dependence of A and the orbital tune ui. Of course, if in addition, A 
depends on extra parameters like the beam energy, then the set Sj is also a function of 
these parameters. Thus if all tori are well-tuned the spin tunes uj s will vary with J and 
any extra parameters. Therefore the distance \u s — m — mi &x — ■ ■ ■ — uj^I from 
the spin-orbit resonance condition ()5.10|) depends on J and any extra parameters. 

Consider a torus for which il(0o) 7^ for each <fr and for each <fr choose a single 
proper UPF [/(•; O ). The corresponding UPRs i/ s (£/(-; O )) define a function of 4> . If 
the torus is well-tuned, then z/ s ([/(-; O )) can be chosen to be a constant. However for 
an ill-tuned torus, v s {U(-; 0o)) usually has to vary with O - This nonconstant function 
can nevertheless be continuous and we present examples in Remarks 13 and 14 in 
Section 6. However, because of the freedom to choose the proper UPRs, z/ s (£/(-; 4>o)), 
this function is usually very irregular and discontinuous. 

If a torus is ill-tuned, then in domains of 4>q where v s {U{-\ 0o)) varies with O > it is likely 
that there are some 4>oS for which the resonance condition u s (U(-; 0o)) = m o + rn ■ u 
holds. 

If the principal solution matrix $ can be written as 



where v{<po) £ [0, 1), -B(0o) is a rea l skew-symmetric matrix with spectrum 
A(£>(</>o)) = {i, —i, 0} and the matrix p(-; (fi ) is in Q(l, lu; d + 1), then p, B, v will be 




(5.11) 
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called generalized Floquet parameters at 0o- In particular f(</>o) will be called a gen- 
eralized Floquet frequency. By definition the generalized Floquet parameter p(9; O ) 
is smooth in 9 and p G 50(3), with p(0; (fi ) = I. Because p(-;0o) has the same 
periodicities as A(-\ 4>q), the p, B, v are indeed generalizations of the Floquet param- 
eters introduced in Section 3. Note that (|5.11jl generates a large class of As, namely 



A{6; <po) = (d$/dff)(6; <j) )$ T {0] <fo), where $ is given by (InTTD . 

If p(9; <f>o), B(4>o), v{4>q) are generalized Floquet parameters at 0o, then from Lemma 
2.1b a W(<po) e 50(3) exists such that W((po)JW T ^ ) = B^), whence $(#;0o) = 
p(0;0o)W / (0o)e Jly(</,o)e (p(O;0o)W / (0o)) T - Thus p(0; 0o)^(0 o ) is a proper UPF at O 
with UPR i/(0 o ). 

Conversely, if [7(0; 0o) is a proper UPF at <po, then by (j2.15j) . $ can be written as in 
(I5.11J1 with p, B, v given by 



Therefore the UPR corresponding to a proper UPF is a generalized Floquet frequency 
and is thus a frequency additional to those in A in analogy with the Floquet fre- 
quency which emerges in the case where A is 27r-periodic. We conclude that the set of 
generalized Floquet frequencies at an arbitrary </> is identical with the set S(0 O ). 

Note that the condition S(0 O ) ^ at some O implies that generalized Floquet param- 
eters exist at that 0o- This can be viewed as a generalized Floquet theorem. Moreover, 
if the torus is well-tuned, then for every <po generalized Floquet parameters exist and 
the set of generalized Floquet frequencies is the same at every <po and is identical with 
the set 5 of spin tunes. 

(9) The spin trajectories at a given O are quasiperiodic if S(0 O ) 7^ but there are A(-] 0o) 
in Q(l,u>; d+ 1) which give rise to nonquasiperiodic spin trajectories (see Section 8). 



6 The Invariant Spin Field and the Quasiperiodicity of 
the Spin Motion 



In this section we present an application of the concepts introduced in Section 5. In particular 
we consider the important situation where for a given A(-, J), a torus admits an invariant 
spin field, a central object in the theory of polarization whose relevance was explained in 
the Introduction. It was pointed out there that the invariant spin field is a solution to the 
T-BMT equation along particle orbits and that it is 27r-periodic in 9 and (f). We now return 
to it and show, among other things, that if the torus is off orbital resonance, the existence 
of a nonunique invariant spin field implies that the system is on spin-orbit resonance. We 
begin by studying general spin fields. 

Let S(9, 4>) be a field such that S(9; O ) := 5(6', O + u9) is the spin of the particle which 
starts with phase 4>q a.t 9 = 0. By (J2.1|) we must have 



p(9;<f> ) := ^;0 o )^ T (O;0o) 
z/(0 o ) := v s {U(-;(f>o)) . 



B(4> ) := t/(O;0 o )^ T (O;0o) , 



d 
d~9 



S(6; O ) = VS(9, uj9 + O ) = A(0, ou9 + 4> )S(9, uj9 + O ) 
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and thus the evolution of the field is described by the partial differential equation 

**k VS = A(9,<p)S . (6.1) 

Here, US (9, 0) = DiS{9, <p)+uj- V^S^O, 0), where D k will denote the derivative with respect 
to the k-th argument, be it scalar or multicomponent. Consider the generalized principal 
solution matrix ip(6 , (ft) , defined by 

V<p = A(9,<f>)<p, <p(0,<f>) = I, (6.2) 

then the general solution of (J6.1)) is given by 

S(6, 0) = f (9, 0)5(0, - oj6) . (6.3) 

This can be seen directly: since T> is a derivation, VS(9,(f>) = ^Vip(9, 0)^5(0, — uj9) + 

y{9,4>)V\S{®,<t>-u)9)^ and£>^S(O,0-cu#)^ = 0. Equations (jHHJ) and (Q are amenable 

to solution by the method of characteristics. In this case the characteristics for (j6.1|) are just 
the spin trajectories S(-,<po) and the characteristics for ()6.2)1 are just the $(-,0o)- 

Now suppose, as in Section 2, that we know one solution of (j6.1|) . Let V(9, 0) be in 5*0(3) 
such that its third column is this solution and make the transformation ip — > ip defined by 

if = Vip. (6.4) 

Then 

Vip = (VV)ip + Wi) = AVip 

so that 

V^j = C v (9, <t>)Ji) , m 0) = V T (0, 0) , (6.5) 
where, in analogy to Section 2, Cy is defined by 

V T (AV - VV) = C V J , (6.6) 

and thus 

C v = -(1/2) Tr[J (V T AV - V T VV)] , (6.7) 
VV = AV - C V VJ (6.8) 

which is analogous to (J2.13|) . 

Equation (|6.5|) also can be solved by the method of characteristics giving 

ip{9, 0) = exp I J [ C v {9', u9' + <j>- cu9)d9' I V T (0, - uB) , (6.9) 



o 
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as is easily checked by direct substitution in (|6.5j) . The generalized principal solution matrix 
thus becomes 



<p(9, 0) = V{9, 0) exp I J I C v {6', lu9' + <P- u9)d9' V T (0, <j> - u9) . (6.10) 



The principal solution matrix of (j5.1|) is easily constructed from the generalized principal 
solution matrix as the following theorem shows. 



Theorem 6.1 Let V(#,0) be a smooth 5*0(3) matrix on a fixed torus, with third column 
satisfying (jfi.lj) . Then the principal solution matrix at 0o on the torus is given by 

**t $(9-0 o ) = V (9,u9 + (f) ) = V(9,l;9 + 4> ) 

■exp I J J C v (9',u9' + <j) )d9'\ V T (O,0 O ) • (6.11) 



Proof: Clearly $(O;0 O ) = I- Moreover, $ satisfies ()5.1j) because (d/d9)V(9,u9 + 



AV - C V VJ and (d/d9) J C v (9',u9' + <p )d9' = C v (9,uo9 + O ). □ 

o 

We now give several definitions. In particular we give a formal definition of the invariant 
spin field. 

Definition 6.2 (ISF,IFF): Consider a fixed torus. 

a) A field S(9, 0) G M 3 is said to be a spin field for ()2.1j) iff S is smooth in (9, 4>), 27r-periodic 
in <p and if it satisfies the partial differential equation (|6.1|) . A spin field S is called an 
invariant spin field (ISF) for ()2.1|) if it is also 27r-periodic in 9 and |<S| = 1. 

b) An SO(3) matrix, V(#,0), is called a frame field for ()2.1|) iff V is smooth in (6,<p), 2tt- 
periodic in and if its third column is a spin field. A frame field is called an invariant frame 
field (IFF) for (|2.1j) iff it is also 27r-periodic in 9. Thus the third column of an IFF is an ISF. 
A uniform IFF, V, is such that the function Cy as defined by ()6.7|) is constant (independent 
of (0,0)) and in [0,1). □ 

Note that in ()6.10|) both V and <p are frame fields. 



Theorem 6.3 a) Consider a fixed torus. If a uniform IFF V exists then, for every 0o, 
?/(■; 0o), defined by U(9; 0o) := V(9, 0o + oj6), is a proper UPF at 0o with v s (U(-\ 0o)) = Cy- 
Moreover the torus is well-tuned and Cy is a spin tune. 

b) Consider a fixed torus. IfV is a frame field then it satisfies the partial differential equation 
(16.81) . where Cy is given by ([6.7)1 . 

Proof of Theorem 6.3a: From Theorem 6.1 it follows that U (•; 0o) is a UPF at 0o with UPR 
v s (U (•; 0o)) = Cy. Of course, U (•; 0o) is in Q(l, u;d+l) so that the UPF U (•; 0o) is proper. 
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Because Cy is contained in each S(0o) we conclude (recall the comment after Definition 5.4) 
that the torus is well-tuned and that Cy is a spin tune. □ 
Proof of Theorem 6.3b: With 

V =: [D 1 , 2 , 3 ] , (6.12) 

and 

W j := (V - A)i) j , (6.13) 
and since V is in SO (3) and S = 3 satisfies (16.1)1 . we have: 

re 1 = (tu 1 • t> 2 )e 2 , w 2 = -(m 1 • o 2 y , ro 3 = , 

so that 

- .4V = (ro 1 • o 2 )VJ r . (6.14) 

Because Tr[J 2 } = -2, it follows from (jgHj) that m 1 • 2 = (1/2) Tr[ J(V T AV - V T VV)}. 

□ 



We now make some further remarks on IFFs and ISFs. 
Remarks: 

(1) The name "invariant frame field" (IFF) is chosen so as to reflect the fact that like the 
ISF it is 27r-periodic in 9. 

(2) If V is an IFF, then we call the o 1 , t> 2 in f!6.12|) pseudo-w 1 , w 2 -axes. They are used in the 
program SPRINT for the numerical calculation of the spin tune El El U21 EH 123 EH| • 
If the IFF is uniform, then o 1 , 2 are called u 1 ,^ 2 — axes (see also [221 121! )• 7 ^- n example 
of a numerical calculation of M 1 ,w 2 -axes can be found in ^3!- See Section 10 also. 

(3) From the proof of Theorem 6.3b it follows for every IFF V that (see also [24 ) 

C v {8, 0) = o 1 ■ (-A(8, 0)o 2 + Vv^j = n ■ o 3 + o 1 ■ (Vv 2 ) . (6.15) 

Thus if O 2 is an ISF, then Cy = and V is a uniform IFF. 

(4) A uniform IFF always exists for the closed orbit so that by Theorem 6.3a the torus 
J = is well-tuned. This follows easily from the fact that in this case proper UPFs 
are uniform IFFs (see Remark 2 in Section 3). Thus at J = every Floquet frequency 
is a spin tune by Remark 1 of Section 5. Note also that for J = any 27r-periodic spin 
trajectory S with \S\ = 1 is an ISF. 

(5) Because the coefficients of the partial differential equation (|6.1|) are 27r-periodic, it is 
plausible that an ISF exists. But as we signaled earlier, it remains as a mathematical 
challenge to prove it and in Remark 14 we will give an example where no ISF exists. 
The coefficients of the partial differential equation (j6.8|) are also 2-7r-periodic so that it 
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is again plausible, but again mathematically challenging to prove, that an IFF exists. 
Note that if an ISF S exists, then an IFF V is easily constructed in analogy with the 
construction of V(9) in Section 3 if there is a 27r-periodic smooth unit vector which is 
nowhere parallel to S (see also |24j). However, if there is no such unit vector then it 
can happen that an IFF V with the third column S does not exist. 

If an ISF does exist, it might not be unique. In that case the torus is on a spin-orbit 
resonance, as shown in the following theorem. 

Theorem 6.4 Consider the case of a fixed torus which is off orbital resonance and assume 
that an ISF, n, exists. Let n be nonunique, i.e. let n be another ISF such that the vector 
product of n and n is nonzero at some (#,0). Then a uniform IFF exists (so that the torus 
is well-tuned) and the system is on a spin-orbit resonance. 

Proof: Let n and n be ISFs whose vector product is nonzero at some (9, 0) and let the torus 
be off orbital resonance. We observe that F := \n x h\ 2 is a smooth and 27r-periodic function 
with T>F = 0. As will be shown below, it follows that off orbital resonance the function F 
is constant. Thus the angle between the ISFs is the same at all points and m, defined by 
m := (n x h)/\n x h\, is an ISF, perpendicular to n. Therefore V := [m x n,m,n] is an IFF. 
Due to Remark 3 the IFF V is uniform with Cy = 0. We conclude from Theorem 6.3a that 
the torus is well-tuned and that is a spin tune. In particular the torus is on a spin-orbit 
resonance (see also Remark 3 in Section 5). 

To complete the proof, we now consider a smooth and 2-7r-periodic function T : M. d+1 — > K. 
for which VT = 0. We define H : R d+1 R via 

H(6, <j>) := F{9, (j) + uj9) . (6.16) 

Hence 7i(8,(j)) is a smooth function 27r-periodic in such that D{K = 0. Because 7i is 
smooth, we have 

H(6, 0) = H(0, 0) + / D 1 H(9 ,(j>)d9 = H(Q,<j)) . (6.17) 
Jo 

Because T is 27r-periodic in 9, we obtain 

= F(2tt, 0) - ^(0, 0) = H(2n, - 2nuo) - H(0, 0) , 

whence by ()6.17j) 

= H(0, - 2tcuj) - H{0, 0) . (6.18) 

Thus for the Fourier coefficients g r := (l/2n) d ■ ■ ■ J Q 27r TC(0, 0) exp(— ir ■ <fi)d(j) of 7~i(0, ■) it 
follows that 

g r = exp(27ur ■ u)g r , (6.19) 
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where r G Z d . Because (1,^) is nonresonant g r vanishes for r 7^ 0. Thus by Lemma 
4.3a H(0,4>) is constant, i.e. independent of 0. Therefore by (j6.17J) H(9,(p) is constant, i.e. 
independent of 9, and then by (j6.16|) T{9, 0) is constant, i.e. independent of 9, 0. □ 

Theorem 6.4 addresses the uniqueness of the ISF as well as its nonuniqueness. In partic- 
ular, the contrapositive of Theorem 6.4 yields: if off orbital resonance a uniform IFF exists 
and if the system is not on spin-orbit resonance, then the ISF is unique up to a sign. This 
behavior was predicted earlier in 



We now focus on the case where an IFF exists and Cy is not constant. We first define 

cy(0;0 o ) :=C v {6,ue + <t>o) . (6.20) 

Since cy(-;0 o ) G Q(l, u; d + 1), its mean, cy(0 o ), an d zero-mean part, cy(0;0 o )> exist. We 
thus have an important decomposition of Cy, namely, 

C v (9, 0) = c v (0 - uo9) + c v (9; (j)-Lu9) . (6.21) 

Since the l.h.s. is 27r-periodic, the r.h.s. is too and in fact cy(0 — uo9) and cy(#; — uj9) are 
individually 2-7r-periodic as we check in Lemma 6.6a below. From (|6.21j) the integral in the 
exponential of (j6.1(JJ) is 

9 

J C v (6', u& + <p- u9)d9' = 9c v ((j) - uj9) + (3 V (9, <p) , (6.22) 



where 





/3 V (0, cf)) ■= J c v {9'- <f> - cu9)d9' . (6.23) 







Clearly Vj3\>(9,(f)) = cy(6 l ; (f> — u9), which leads to consideration of the partial differential 
equation 

Va(9, 0) = c v {9- cj) - u9) =: f(9, 0). (6.24) 
Then for every solution a of (J6.24j) 

(3 v (9,(f)) = a{9,(t))-a{0,<j)-u9) . (6.25) 

The existence of a 27r-periodic a will be important below (note that T is always 27r-per iodic). 
We now write 

cy(0o) = M0o) + k v ((po) , (6.26) 

where the integer fcy(0 o ) is uniquely determined by the condition z/y(0 o ) G [0, 1). The 
generalized principal solution matrix from (|6.1U|) now becomes 

= V(0, <p) exp (j{a{9, <p) - a(0, <p - uj9) + k v (<p - uj9)9 + vy{4> - lu9)9)^J V t (0, <p - uj9) 
= U(6, 0) exp ( JV v (0 - u9)9) W T (0, <f>-u6), (6.27) 
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where 



U(9, 0) := V(0, 0) exp ij(a(9, 0) + fc v (0 - "6)9) J . (6.28) 

Then the principal solution matrix becomes 

$(0; O ) = U(6; O ) exp (^(00)0) U T (0; O ), (6.29) 

where 

[7(0; 0o) :=U{6,u;6 + <f>a) 

= V(0, cu0 + O ) exp ^(a(0, cu0 + O ) + k v (<p )9^J . (6.30) 

We can now state and prove the next basic result of this section. The proof will depend, in 

part, on Lemma 6.6b which follows later in order not to break the flow. Theorems 6.5c-d 
use the Diophantine condition. 



Theorem 6.5 Consider a fixed torus Jq. 

a) If an IFF, V, exists and (I6.24|) has a smooth and lis -periodic solution a, then [7(-;0o) 7 
defined by fl6.3()j) . is a proper UPF at O with UPR equal to z/y(0 o ). 

b) If the conditions of Theorem 6.5a hold and if the torus is off orbital resonance, then U, 
defined by (I6.28J) is a uniform IFF and Cy = vv- Moreover the torus is well-tuned and 
[/(■;0o), defined for every O by ()6.30p . is a proper UPF at O whose UPR is a spin tune 
and u s (U(-; O )) = C u = u v . 

c) Let < t < r — d — 3 and let A(-, Jo) be in C r . If an IFF, V , exists in C r and if 
uj(Jq) G ^(t) then a uniform IFF exists (and thus the torus is well-tuned). 

d) Let r > 2d + 4 and let A(-, Jo) be in C r . If an IFF, V, exists in C r for every u(J ) 
in a Borel subset 1Z of Mr then a uniform IFF exists (and thus the torus is well-tuned) for 
[i-almost every to{ Jo) in 71. 

Proof of Theorem 6. 5a: It is clear from ()6.30|) that U (•; 0o) is in Q(lj uj;d+l) and the result 
follows from fl6~2"9"|) . □ 
Proof of Theorem 6.5b: From (J6.28)) . U and V have the same third column and U is smooth 
and 27r-periodic. Therefore U is an IFF. From (J6.7j) an easy calculation gives Cy = z/y. Thus 
by Lemma 6.6b, Cy is constant, independent of 0o- Hence U is a uniform IFF. Therefore 
by Theorem 6.3a [/(•; 0o) is a proper UPF at O with UPR z/y which is a spin tune and the 
torus is well-tuned. □ 
Proof of Theorem 6.5c: From ()6.7)1 . Cy is in C r ~ x since V and A are in C r . Using the fact 
that the torus is off orbital resonance (since uj(J) G Q(t)), we have by Lemma 6.6b that 
cy is a constant. Thus, by ()6.21|) . T = Cy — cy is in C r_1 . It follows from the condition 
r < r — d — 3 and Lemma 4.7 that 

a(9,(f))= — -J" m exp(im • (0,0)) , (6.31) 

^— ' im ■ (I.lo) 

m£j, d + 1 \{0} V ' ' 
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is smooth and 27r-periodic in 9 and <fi and that it satisfies (|6.24|) . where f m denotes the m-th 
Fourier coefficient of T . The claims now follow from Theorems 6.5b and 6.3a. □ 
Proof of Theorem 6.5d: The interval (d+ 1, r — d — 3) is not empty so that we pick a r in 
that interval. Because r < r — d — 3 we have, by Theorem 6.5c, for u{ Jo) £ (fi(r) n TV) a 
uniform IFF (and thus a well-tuned torus). Because r > d + 1 we have by Lemma 4.8 that 
/i-almost every tu( Jo) in 72 is in fi(r) n 72.. This proves our claim. □ 

We now complete the discussion by stating and proving the Lemmas 6.6a and 6.6b mentioned 
earlier in this section. 



Lemma 6.6 a) IfV denotes an IFF, then cy(0 — oj6) and cy(0; (p — oj9) are lit -periodic in 
9 and in <fi. 

b) If V denotes an IFF and if the torus is off orbital resonance, then cy,i/y and fcy are 
constant. 

Proof of Temma 6. 6a: The periodicities in ()6.21|) can be demonstrated as follows: 
c v (0 - u(6 + 2tt)) = lim 1 / C v (9', <p - uo{9 + 2tt) + uo9')d9' 







= lim - / C v (9' - 2ir, <f> - uo(9 + 2tt) + uo9')d9' 
= lim - / C v {9",(f)-u9 + u9")d9" 

T^oo T J_ 2n 

= c v (0 - w0) + lim i f / - / )c v (9", <p-u9 + u9")d9" 

T^oo 1 VJ-27T JT-2-kJ 

= c v {(j) - uj9) , (6.32) 

where at the second equality we have used the fact that Cy(-,0) is 27r-periodic and at the 
last equality we have used the fact that Cy is bounded. This shows that the first term on the 
r.h.s. of (j6.21|) is 27r-periodic in 9 and thus that all three terms in ()6.21|) have this periodicity 
property. That all three terms in (j6.21J) are 27r-periodic in is trivial. □ 
Proof of Temma 6.6b: Off orbital resonance we find, due to ()6.20|) and by applying Lemma 
4.3c, cy(0) = cy(0 o ) = (l/27r)( d+1 ) J^---J^ C v d9d<\>. Thus c v , v v and k v in (I6~2T)1) are 
(/^-independent in this case. □ 

Remarks: 

(6) The 27r-periodicity in 9 and of cy(<f) — uo9) and cy(6 l ; tp — uo9) is suggested by examining 
the formal Fourier series 

E c m ,„e i(n0+m ^ ) of C v where neZ,meZ d . Then it is easy to 

m,n 

show that the resonant module part of this sum, defined by m- (1, u) — for m £ Z d+1 
corresponds to cy(0 — uj9). The remaining part with m ■ (1,uj) ^ corresponds to 
cy(^;0 — uj9). Their formal Fourier series display the 27r-periodicity in 9 and 0. Off 
orbital resonance the relation n+m-u = implies that n = m = so that cy(0o) — co,q, 
i.e. it is independent of 0o as in Lemma 6.6b. 
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(7) Under the conditions of Theorem 6.5a, z/y(0 o ) is the UPR at 0o associated with the 
UPF defined in (|6.3U|) . However the UPF is not unique since the principal solution 
matrix in ()6.29j) can be written as 

$(0; O ) = U(6; O ) expWfo)) exp(JMh)O) exp(-Jh(<j)o))U T '(0; O ), 

so that £/(#; 0o) = 0o) exp(j7/i(0o)) is also a UPF for an arbitrary smooth h(4>o). 

(8) Under the conditions of Theorem 6.5a, and using the notation in (j6.12J) . S(9;(p Q ) = 

O )D 3 (O, 0o ) gives a solution of (j2.1j) which is in Q(l,u;d + 1). This is easy to 
check by (|6.29jl : since V T (0, 0o)t> 3 (O, 0o) = (0,0,1) is the eigenvector of J with zero 
eigenvalue, S(9;4>o) — ^> 3 (9, (f)(9)) an d ^vi^o) has dropped out. By ()6.29j) all other 
linearly independent solutions are in Q(l, lu, z/y(0 o ); d + 2). 

(9) The IFF underlying the definition (j6.26j) of z>v(0o) is, of course, not unique. For 
example, if z/y(0 o ) G (|, 1) then by changing the signs of t> 2 and 3 we find Vy{4>o) G 
(0, 5). Thus in analogy to the case J = in Remark 7 in Section 3 we can choose V 
such that z/y G [0, |]. 

(10) If V denotes a uniform IFF, then cy, cy, z/y, &v are constant and cy is zero. In particular 

z^y = Cy. 

(11) Under the conditions of Theorem 6.5b it is possible, by rotating an arbitrary IFF into 
a uniform IFF, to construct a proper UPF at every point 4> of the torus such that the 
UPR is independent of (fi and is a spin tune. 

(12) Theorem 6.5a holds on and off orbital resonance although if uj does not satisfy a 
Diophantine condition, appropriate solutions of (j6.24j) may not exist. Note also that 
in Section 8 we will consider an example where the conditions of Theorem 6.1 hold 
but where it turns out that the conditions of Theorem 6.5a cannot hold due to the 
presence of nonquasiperiodic spin motion. 

(13) To illustrate Theorem 6.5, it is instructive to consider simple, but perhaps unphysical, 
models. Here we will use the model defined by A(6, 0) := y/2J cos(0 — 6) J with d = 1. 
This represents a precession around the vertical and, of course, such an A will not 
be found in real storage rings. Note that A is a smooth and 27r-periodic function 
and, because (0, 0, 1) is an ISF, we can choose the IFF V as V — I so that by (J6.7J) 
Cy(#, 0) = V2J cos(0 — 9). We first consider a case on orbital resonance, uo = 1. From 
(16.20(1 . cy(6*;0 o ) = V2J cos(0o), whence cy(0o) = cy(#;0o), cy = 0. Clearly ()6.24() is 
solved with a = and ()6.30(] gives U(9; 0o) = exp (J7"/Cy(0o)6*). This [/(•; 0o) is a proper 
UPF at 0o with UPR z/y(0 o ), where z/y(0 o ) and the integer fcy(0o) are determined by 
the condition 

&v(0o) + ^v(0o) = V2J cos(0 o ) . (6.33) 

Obviously, if J > 0, E has uncountably many elements so that the torus is ill-tuned; 
in particular, by Theorem 6.3a, no uniform IFF exists. 
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We now consider the case of irrational uj. From ()6.20|) . cy(0; 0o) = V%J cos((u; — 1)0 + 
0o), thus Cy = and Cy(0; <f> ) = Cy(0; (fi Q ). Since c v = 0, v v and ky are zero and since 
Va(6,(f)) = c v (0;0 - u9), we have a(9,<p) = V2J sm(<p - 6)/(u - 1). Then 

gives U(9, (p) = exp ^J(y/2J / (uj — 1)) sin(0 — 0)J . One observes that U is a uniform 

IFF and = 0. Hence by Theorem 6.3a the torus is well-tuned and on a spin-orbit 
resonance. 

(14) In Theorems 6.3-5, we assumed that an ISF exists and for the physically interesting 
.As we hope this is true. It is also clear from Section 3 that for J = an ISF does 
exist. Nevertheless, for J / this assumption in the theorems is not superfluous as 
we will now show by constructing, for a fixed torus, an A such that proper UPFs exist 
at every O but an ISF does not. 

We consider the case where d — 1, ui — 1 and where A(6, 0) depends only on <p — and 
is given by 

/ o -n a ((f>-9) n 2 (0-0) 
A(o, (fy = A{(f> - o) := n 3 ((f>-e) o -n x {<t>-&) 

\-O 2 (</>-0) fii(0-0) o 

Here the function Q is smooth and 27r-periodic and we assume \Q\ < 1. The principal 
solution matrix is given by $(0; <fio) = exp(^4(0 o )0). Because A is skew-symmetric 
with spectrum \(A((fio)) = O'|^(0o)|> — i | O(0 O ) 1 , 0}, it follows by Lemma 2.1b that a 
50(3)-matrix W(cf>o) exists such that 

A(<Po) = iy(0o)|^(0o)|^ T (0o) • (6.34) 

Hence 

$(0;0o) = iy(0o)exp( l 7|fi(0 o )|0)H/ T (0 o ) , (6.35) 
so that W((f)o) is a proper UPF at <p Q with UPR |f2(</> )|. 

We now proceed to construct an Q such that an ISF does not exist. If an ISF n exists 
then, by ()6.3|) and Theorem 6.1, we have 

$(2tt; (f> )n(0, O ) = n(0, O + 2tvuj) , 

for all 4> . Then, by (|6.35j1. 

exp(J\n( ( f ) o)\2-K)W T ( ( j ) o)n(0; ( j ) o) = W T (<J) )n(0; O ) • 

If |f2(0 o )| 7^ integer, i.e. if f2(0o) 7^ 0, then every eigenvector of exp(j7"|fi(0o)|27r) for 
eigenvalue 1 is parallel to (0,0,1). In addition, since A(4>o)Q(4>o) = 0, ()6.34j) implies 
that 

jw T ((f> )n(<j> ) = o . 
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It follows that W T (4> )n(0; O ) and W T (4> )Q(4> ) are parallel to (0, 0, 1). Hence n(0; O ) = 
const x fl((f>o). We conclude that for fi(0o) 7^ 



n(O;0o) = ±n(0o)/Wo)|- 

Thus if an ISF n exists, then 



n(0, 0o 



±fi(0o)/|^(0o)| if ^(0o)^O 
f if fi(0 o ) = , 



where £ is arbitrary of norm one. To obtain our example, we present a smooth Q such 
that fl(<po)/\fl(<p )\ does not have a continuous extension to those <f) , where fl(4>o) = 0. 
Let ip(x) := 1/ sin(x) and 

OU ^ - / (l/2)sin 3 (0o)(cos(^(0o)),-sin(^(0 o )),O) if sin(0 o ) ^ 
"^-\0 if sin(0 o ) = O. 

It is easy to show that Q is smooth with < 1 so that A(6,<p) is smooth and 27r- 
periodic. However 

fi(0 o )/|fi(0 o )| = | Sm( .^ | (cos(^(0o)),-sin(^(0 o )),O) , sin(0 o ) ^ 
I sm(0o)| 

does not have a continuous extension to all of 0o- Thus there is no ISF. 

In summary, we have an example of a torus on orbital resonance for which a proper 
UPF exists at each O but which is ill-tuned because the set {(1/2)| sin 3 (0 o )| : O e M} 
is a subset of H and has uncountably many elements. Furthermore there is no ISF and 
thus no IFF. Of course, an A as exotic as that defined here will not emerge from the 
fields of a real storage ring. 

(15) Now that IFFs are available we present another generalization of the Floquet Theorem. 
Starting with ()6.29|) and assuming that the conditions of Theorem 6.5a hold, we write 

$(0; O ) = U{6- O ) exp(^v(0o)fl)t/ T (O; 0o) 

= U(6; <P )U T (0; <f> o )U{0; O ) exp(^ v (0 o )^)f/ T (O; O ) 

= p(0; O ) exp (5(0 o >(0 o )0j , (6.36) 

where 

p(0; O ) := 17(0; 0o)f/ T (O; O ) , Wo) := ?7(0; O )^ T (O; O ) , K0o) := M0o) , 

and where U is given by (|6.30j) . Thus U is a proper UPF. It follows from Remark 8 
of Section 5 that p(-; 4>o), B((j) ), v(<fio) are generalized Floquet parameters at <fi . Thus 
under the conditions of Theorem 6.5a we have generalized Floquet parameters and 
have obtained a generalization of the Floquet theorem at every <f) . Note that ^y(0o) 
is given by (|6.26j) . 
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(16) For speed and practicality, simulations of particle motion in storage rings are often 
made with the approximations that the fields at the ends of magnets fall to zero 
abruptly (the "hard edge approximation") or that a magnet has zero length but the 
correct field integral (the "thin lens approximation"). Then the .As are not smooth in 
9 so that the theory of this paper would have to modified. For example an ISF would 
be defined without imposing smoothness in 9. The IFF would also not be smooth in 
9. We return to these matters in Section 10. 

We now prove a simple partial converse of Theorem 6.3a which gives a large class of well- 
tuned .As with a uniform IFF. 

Theorem 6.7 Let q(9, <fi) be a SO(3) matrix which is of class C 2 and 2n -periodic such that 
q(0, 4>) = I and let b be a constant real skew-symmetric matrix. Let A be defined by (|5.2j) 
with A = (Vq + qb)q T . Then a uniform IFF exists, so that the torus is well-tuned. 

Proof: Because b is skew-symmetric, by Lemma 2.1b a constant 5*0(3) matrix W exists 
such that b — {y + m)WJ"W T , where v e [0, 1) and the integer m is constant. Defining 
V := qW exp(j7"m#), we observe that the smooth 5*0(3) function V(9,<f)) is 27r-periodic and 
that T>V = AV — vVJ. Hence Cy = v. For the third column of V we have 

£>t, 3 = £>V(0,0, 1) = AV(0,0,1) . (6.37) 

So it is an ISF and therefore V is an IFF. The IFF V is uniform so that by Theorem 6.3a 
the torus is well-tuned. □ 

Remark: 

(18) We obtain the principal solution matrix as in (J6.29)) . with z/y = v and 

U(9; (f> ) = q(G, u9 + O ) W exp (Jm9) . (6.38) 
Combining (j6.29j) and ()6.38|) gives the generalized Floquet form 

$(#; O ) = q (9, u9 + O ) W exp (J(m + v) 9) W T = q{9, u9 + O ) exp(W). 

We now apply some ideas of this section to some more simple models. Then in Section 
9 we return to the general case and consider the spectrum of the principal solution matrix. 

7 The Single Resonance Model 

In this section we examine the so-called "single resonance model". This model provides 
an approximation for the matrix A and it is frequently used in approximate descriptions 
of spin-orbit resonance effects both in proton and electron rings [HIEE!- The reason for its 
popularity is that it can be solved exactly and delivers useful indications about behavior 
near spin-orbit resonance even though A is an approximation. Since it can be solved exactly 
it provides an example for O ) m Q6.29J) and (|6.30p . 
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The single resonance model effectively describes the spin motion for a particle limited 
to a harmonic vertical betatron oscillation around a horizontal circular design orbit. Thus 
there is only one pair of action-angle variables J, <fi, and d — 1. The matrix A then contains 
terms due to a constant vertical magnetic field and terms due to a radial quadrupole field 
oscillating with the tune to = d<p/d9. The radial field can be decomposed into two counter 
rotating horizontal components rotating with the tune ui. Close to resonance it is sufficient 
to consider just one of the rotating field components and to neglect the other ^SJ HH]- The 
matrix A then takes the form: 



-k-kk A(4>) \- 



/ —o"i cr 2 v2J sin(0) 

o-i -a 2 v / 2Jcos(0) | . (7.1) 

\ —CT2V2J sin(0) 02V2J cos(0) 



where <j\ and <Ji are real constants describing the strength of the vertical field and of the 
rotating horizontal field respectively. Note that A is a smooth function of 0, independent of 
6. 

Our first aim is to show that a proper UPF exists at every O - We first consider the case 
a 1 7^ u), (T2V2J 7^ and define a := a/ {a\ — uj) 2 + 2o\ J. Then it is easy to verify that (see 
also 0) 



1 / CX2V 2Jcos(</>) 
O 3 (0) = - <7 2 y/2Jsm(<f>) I • (7.2) 



o x — UJ 



This is 27r-periodic and it satisfies (j6.1|) . It is thus an ISF. Since (1,0,0) x O 3 (0) 7^ 0, we 
can define t> 2 (0)_Lo 3 (0) by 



2/ ^_ (1,0,0) xo 3 (0) 



» <*> = = v^mm [ 1 ' (7 ' 3) 



and O 1 (0) by D 1 := 2 x o 3 . Thus V := [o 1 , 2 , t) 3 ] is smooth and 27r-periodic. It is thus an 
IFF. By (IFrToTl .(l7^1). (l73|) we obtain 

Cy(0) = m • O 3 (0) + ^(0) • (lV(0)) = a + g2 ^Jcol^) ' (7 - 4) 



so that 



c2-k 

c v (0o) = hm (1/T) / d6' C v (uj9' + (j) ) = (l/27t) I dcj> C v (<p) , (7.5) 



where we use the 27r-periodicity of Cy. Thus cy(0 o ) is independent of O - Since jSUl p. 379] 



71 



(oi — u)) 2 sin ((f)) + a 2 cos 2 (0) 2cr|(Xi — k>| 
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and since by ()7.4|) and (|7.5|) 

o<jj(o\ — to) '" 2 ' 

c v = a + 



2n J (a i — uj) 2 cos 2 (0) + a 2 sin (0) 
auj(<Ji — uj) 



71 Jo (ci — c<j) 2 cos 2 (0) + cr 2 sin 2 (0) 
2acj((7i - u) rw/2 

71 Jo — sin 2 (0) + cr 2 cos 2 (0) 



Z / , _9_„9/^ » ( 7 - 6 ) 



we can make the assignments 

c v = a + w sign((Ti - u) , c v (</>) = a 2 _ 2o - 2 J cos 2 (0) ~~ ^ ^S 11 ! "! ~ w ) • ('•') 

Because cy does not depend on 9 one finds that ()6.24|) has smooth solutions a((f>) which do 
not depend on 9 and that one of those solutions reads as 

a((f>) = - [ #'cv(0') . (7.8) 



Since = J Q 27r d<pcy((j)) and by the 27r-periodicity of cy(0) it follows by (|7.8|) that a(0) is 
27r-periodic in 0. Thus Theorem 6.5a applies and one concludes that [/(•; 0o), defined by 

d0'cv(0') + M) , (7-9) 



U(9; O ) = V(0 O + uj9) exp / ' 



is a proper UPF at 0o with UPR z/y. Here z/y and the integer ky are uniquely determined 
by cy via ()6.26j) and (|7.7j) . Because z/y is contained in each S(0 O ) we conclude (recall the 
comment after Definition 5.4) that the tori are well-tuned when o\ ^ uo and o^v^J 7^ 0. 
We now discuss the general case using the machinery of Section 6. 

Proposition 7.1 The single resonance model has a uniform IFF for every value J of the 
orbital action variable. Hence the corresponding torus is well-tuned. 

Proof: For y{9,<f>) := exp(0j7") one obtains 

y T {Ay -vy) = E , (7.10) 

where 

/ -ax+u 
E := 0-1 - u -0-2V2J 
\ ^v 7 ^/ 

This has eigenvalues X(E) = {io~, —ia, 0}. If a = 0, then £7 = and, due to ()6.6|) and ()7.10|) . 
y is a uniform IFF and Cy = 0. Thus by Theorem 6.3a the tori are well-tuned and is a 
spin tune. If a ^ 0, then by Lemma 2.1b we can choose W G SO (3) such that 
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where the integer k is chosen such that v G [0, 1). We now define 

U(6;<P) := y(9;4>)Wexp(kej) , 

and obtain 

U T {AU-VU) = vj . 

It follows by ()6.6|) that U is a uniform IFF and Cy = v = a — k. Thus by Theorem 6.3a the 
tori are well-tuned and a — k is a spin tune. □ 

Remarks: 

(1) The proof of Proposition 7.1 and Definition 5.5 show that the spin tunes uj s of the 
single resonance model have the form 

uj s = ecr + ju + k , (7-11) 

where e G {— 1,1}; j, k G r L. Conversely every constant of the form (|7.11|) is a spin 
tune of the single resonance model, if it is in [0,1). In particular the set S(0o) is 
independent of 0o- Note that for the single resonance model, spin tunes exist also on 
orbital resonance, i.e. for rational u. See Remark 4 too. 

(2) The case oi\[2J = represents the absence of betatron motion, i.e. motion on the 
design orbit. In this case s,j, k can be chosen in (j7.11j) such that the spin tune reduces 

to 0\. 

(3) From the expression for a it is clear that during variation of a\, the spin tune in (|7.11j) 
comes closest to the spin-orbit resonance uo s = k + jut when o\ = u. However for the 
case o"2\/2J ^ of principle interest, the resonance condition is not reached. 

(4) If u} is an integer, the matrix A is one-turn periodic and the ISF can be obtained 
as the eigenvector of length 1 of the one-turn principal solution matrix. Moreover, 
every proper UPF is one-turn periodic and its UPR can be obtained from the complex 
eigenvalues of the one-turn principal solution matrix just as in Section 3. Of course, this 
UPR is a spin tune. If u is rational, the ISF can be obtained as the eigenvector of length 
1 of the appropriate multi-turn principal solution matrix. See, e.g. |51j . Every proper 
UPF is then multi-turn periodic and its UPR is extracted from the corresponding 
complex eigenvalues and is again a spin tune. Note that this circumstance that a spin 
tune exists even on orbital resonance has its origin in the facts that the single resonance 
model has only one orbital frequency and that A is independent of 9. 

(5) If the orbital tune is rational then with Proposition 7.1 we have an example of a torus 
which is on orbital resonance but is nevertheless well-tuned. This is an example for 
which the torus may be on orbital resonance but still satisfy the conditions of Theorem 
6.3a. 
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8 The Moser-Siegel Model 



In this section we construct and study an illustrative but unphysical model which, for reasons 
which will become clear, we will call the "Moser-Siegel model". This model has d — 1, can 
be solved exactly, and has two real parameters ai,a 2 , where o\ > 1. For certain choices of 
the orbital tune u and of a±, a 2 the Moser-Siegel model provides an example of spin motion 
which is nonquasiperiodic. The matrix A of the Moser-Siegel model takes the form: 



**k A:=A 2 iJ, A 2 i(0,<f>) ■= Yl a i k l Vkisin(k(f)-W) , (8.1) 

where 



k,l=l 



1 , if ka 2 - I > 

-1 , if ka 2 — I < 



and by definition A is smooth (in fact C°°) and 27r-periodic. The latter follows from the 
convergence of X/feLi k n <Ji k for a nonnegative integer n, which follows from the ratio test and 
which implies that the series in (J8.1)) and the series of all its derivatives converge uniformly 
(see Hi 8.6.3], jUl p.117]). 
For 0o = 0, 

00 

A(9; 0) = f(9)J , f(6) = £ a^ k - l Vkl sm((ku - 1)6) . (8.2) 

k,l=l 

Clearly / G Q(l,u; 2) and for an irrational u, f = by Lemma 4.3c, so that 

mO)=exp(Jg(9)) , (8.3) 

where 

e 

9(0) := J ROW ■ (8-4) 


Because the function g has been used in |52[ Paragraph 36], we call our model the Moser- 
Siegel model. 

Proposition 8.1 For some uj ^ Q there exist values ofa\, a 2 such that g is unbounded. For 
these values, exp(ig) is not quasiperiodic whence $(-;0) is not quasiperiodic. 

Proof: It is shown in Paragraph 36] that values of o\ , o~ 2 and of to Q exist such that g 
is unbounded. Thus, for these values, / is an almost periodic function whose integral, g, is 
unbounded. It then follows (see [HI Chapter 6]) that exp(ig) is not almost periodic, whence 
at least one of cos(g) and sin(g) is not almost periodic. □ 

Remarks: 

(1) Proposition 8.1 shows that for certain values of <7i, a 2 and of u the principal solution 
matrix at 0o = is not quasiperiodic, so that H(0o = 0) = 0. In particular, for those 
values the torus is ill-tuned so that, by Theorem 6.3a, it has no uniform IFF. 
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(2) The unit matrix obviously provides an IFF for the Moser-Siegel model. Thus Propo- 
sition 8.1 demonstrates that the existence of an IFF is neither sufficient for having a 
uniform IFF nor for having a well-tuned torus. In the context of Theorem 6.5a this 
means that the existence of an IFF does not necessarily admit a smooth solution a of 
()6.24j) . which is 27r-periodic. 



9 Spectral Analysis of Quasiperiodic Spin Motion 

In the code SPRINT, ISFs are calculated nonperturbatively in three ways, namely by stro- 
boscope averaging, by the SODOM2 algorithm and by adiabatic antidamping [HI El El E21 
[EE I2H1 122] ■ Perturbative algorithms for obtaining ISFs are listed in the Introduction. More- 
over the spin tune is calculated nonperturbatively in SPRINT by logging the spin precession 
around the ISF or by the SODOM2 algorithm. We will make further comments on these 
simulations in Section 10. 

The fact that under appropriate conditions the set of generalized Floquet frequencies is 
the set of spin tunes, suggests a further way to obtain the spin tune, namely by spectral 
analysis of quasiperiodic functions. The way forward is contained in the following theorem. 
As we shall see, we can also use spectral analysis to construct the ISF. 



Theorem 9.1 Consider a uniform IFF V on a fixed torus. Then the following holds. 

a) Let S be a spin trajectory at 0o- //0 6 [Cy], where the equivalence class [Cy] is defined in 
Section 5, then S is in Q(l,u;d+ 1). 

b) Let S be a spin trajectory at O a> n d let the torus be off orbital resonance. Also, let (jL [Cy] . 
Then S, defined by 

S(9) := lim N Awm exp(i6*m • (1, uj))ai S, m ■ (1, uS) I , (9.1) 

iV— >oo ' J \ I 

\\m\\<N 

with AN, m '■= rin=i N+ n+T"^ ' sa ^ s fi es ^ e relation 

S{6) = n(9, 0o + w6)S(Q) ■ n(0, O ) , (9.2) 

where n is the ISF n := V(0, 0, 1). Moreover, S is a spin trajectory in Q(l, u; d + 1). 

c) For arbitrary <fio, 

★*t A($(-;0 O )) C {eCy + m- : (e,m) £ {0,1,-1} x Z d+1 }. 

Proof of Theorem 9.1a: For fixed 0o> define U(9) := V(9, 0o + w $) so that by Theorem 6.3a 
U is a proper UPF at 0o an d we have 

$(0; O ) = U{9) exp(JC v 9)U T (0) . (9.3) 

Thus if G [Cy], then $(•; (fi ), and therefore S, is in Q(l, u;d+l). □ 
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Proof of Theorem 9.1b: Because the U, denned in the proof of Theorem 9.1a, is in Q(l, u; d+ 
1) and because the torus is off orbital resonance we can use Lemma 4.3b-c to write 

U(9) = Hm V A N>m exp(i6m ■ (l,u))a[U,m- (l,u) ) . (9.4) 

N— >oo ' J \ I 



||m||<iV 
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A± := \ | =B 1 | , A := | | , (9.5) 

and recalling (j2.12j) . we have 

exp{JC v 6) = A+ exp{iC v 6) + A exp(-iCytf) + A , (9.6) 

so that 



ayU exp(j7"Cv -); m " (Ij^J = a J ^UA + ,m ■ — Cy j + ayUA_,m ■ (1,uj) + Cy 

+a^C/A ,m- (l,w)^ . (9.7) 

Since ^ [Cy], then ±Cy ^ {m • (l,u)) : m G Z d+1 }. It follows from Lemma 4.3e that 
±C V £ A(t7). Then = a\U,m ■ (1, w) =F Cv^) so that (jHHJ) leads to 

aHJexp^Cy -),m- (l ( w)j = a\UA , m ■ (1, cuU . (9.8) 
Combining (JHEJ), (jOjl . we obtain 

U(8)A = lim N Azvm exp(i0m ■ (1, cj))a( C/A , m ■ (1, u;) 

\\m.\\<N 

lim y A/v,m exp(i6*m • (1, a;))a( t/exp( l 7Cy •), m • (1, a 



lim N Awm exp(z0m • (1, oj))a[ <&(•; </> )£7(0), m • (1, u) 



Hence 



N^oo 

\\m\\<N 



U(9)A U T (0)= lim V 4, m exp(fc-(l, W ))a($(^o),m-(l,w)). (9.9) 

m gZii+l 
||m||<JV 
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Then, because 5 is a spin trajectory at 4>o, 

U(9)A U T (0)S(0) = lim V A Nm exp(i0m • (1, w))o( $(•; 0o)5(O), m ■ (1, w) 



JV^oo 

||m||<AT 



lim N AN m exp(i9m ■ (l,u))a[ S,m ■ ) 

m«=y.<i+l v 7 



||m||<JV 



5(0) , (9.10) 



where the last equality follows from flHIEJ). From (jZHj), O, fl2ZEDl) it follows that 

5(0) = f>(0)A o f/ T (O)5(O) = (a(0;0 o )[/(0) -C V U{Q)3^ A f/ T (0)5(0) 

= A(0; o M0)A o t/ T (O)5(O) = O )5(0) . (9.11) 

Thus 5 is a spin trajectory at 0o- By (|9.1()j) and because U is proper, 5 is in Q(l, u; d + 1). 
If 5 is in Q(l, a;; d + 1), then by Lemmas 4.3b and 4.3c 

5(0) = lim V Aiv m exp(^m- (1, u;))a( 5, m • (l,w) J , (9.12) 

||m||<JV 

so that in this special situation, 5 = 5. Thus, for an arbitrary spin trajectory 5 at (fio, the 

double transform 5 of 5 from the double application of ()9.1j) is equal to the single transform 
5 so that (jOTIjl yields 

£/(0)A o f/ T (O)5(O) = 5(0) = 5(0) . (9.13) 

If0 = 0, (jniSl) becomes an eigenproblem for 5(0). The solution is 5(0) = ±|5(0)|£/(0)(0, 0, 1). 
Inserting this into ()9.13|) yields 

5(0) = ±|5(O)|n(0, 0o + u6) = n{9, 4> + u9)S(0) ■ n(0, O ) , (9.14) 

where n denotes the third column of V and where 5 is an arbitrary spin trajectory at <J)q. □ 
Proof of Theorem 9.1c: By f!9.3j) and (|9.fij) we have 

$(0 ; O ) = £/(0) exp(iC v 0) + A_ exp(-zC v 0) + A ^ U T (0) . (9.15) 
If A G A[UA + exp(iC v -)U T (0) j , then 



A - C v G A f/A+f/ T (0) . (9.16) 



Because £/ is proper and with Lemma 4.3e, A(UA + U T (0)) C {m ■ (1, o>) : m G Z d+1 }. Thus 
with (fnHEJ), A G {C v + m • (1, u) : m G Z d+1 }. Then 

A^A + exp(iC v -)^ T (0)^ c{C v + m-(l,w):mGZ w }. (9.17) 
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In an analogous way 

A^f/A_exp(-iCv)^ T (0)^ c{-C v + m-(l,Lu):meZ d+1 }. (9.18) 

Moreover, 

A^C/A tf T (0)^ C {m ■ (l,w) : m G Z d+1 } . (9.19) 

Combining f|9~T5j) . (j9~T7j) . (EHHJ) . (EH3j) gives 
A($(-;0 O )) 

C ^A^A + exp(zC v -)^ T (0)^ U A A_ exp(— zCy -)?7 T (0)^ UA^f/A o [/ T (0) 

C {^C v + m • (1, u) : (e, m) G {0, 1, -1} x Z d+1 }. □ 
Remarks: 

(1) Let the conditions of Theorem 9.1b hold and let the spin trajectory S at 0o be in 
Q(l,u;d+ 1). In this special situation ()9.12|) holds, i.e. (|9.1j) becomes the spectral 
expansion of S. However, in general S is not in Q(l,u; d + 1), i.e. in general (J9.1)) is 
not the spectral expansion of S because only the tunes m • (1, uS) appear in ()9.1|) . 

If the conditions of Theorem 9.1b hold then S is parallel to an ISF and with (|9.1|) and 
()9.2j) one could, at least in principle, attempt to compute the ISF by doing numerical 
spectral analysis on an arbitrary spin trajectory. Of course, by Theorem 6.4, the ISF 
n is unique up to a sign and, by Theorem 6.3a, the torus is well-tuned. 

(2) Consider a fixed torus and assume that a uniform IFF exists so that, due to Theorem 
6.3a, the torus is well-tuned. Then, due to Remark 8 in Section 5, the set of generalized 
Floquet frequencies is the same at every O an d is identical with the set H of spin tunes. 
Therefore Theorem 9.1c implies, for arbitrary O > that for every A in A($(-; 0o)) \ { m ' 
(l,o;) : m G Z d+1 } the fractional part of A is a generalized Floquet frequency at O > 
and in particular a spin tune. Thus as conjectured earlier, spin tunes can indeed be 
obtained by spectral analysis. 

Theorem 9.1c addresses the spectrum, A($(-; O ))> f° r the conditions stated but gives 
no information on its dependence on <f) . However, under certain conditions, the special 
parameter dependence of A given in (|5.4j) guarantees that A($(-;0 O )) is independent of 4> 
as we show in the next theorem. 

Consider a fixed torus and denote $(0 + 27rA^; <fi )§ T (2tiN; o ) by H(0; <f> Q ), where is 
an integer. We conclude from (|5.1|) and (|5.4|) that n satisfies the initial value problem 

= A(fl;27riVW0 o )n(Mo), n(O;0 o ) = J. (9.20) 

Because $(0; 2'kNuj + <f> ) is the unique solution of this initial value problem (see also Section 
2), we obtain 

$(9; 2itNuj + O ) = $(0 + 2ttA^; o )® 7 \2ttN; o ) , (9.21) 
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valid for arbitrary 0o and arbitrary integer N. Thus the basic property of A in ()5.4|) leads 
to the basic property of $ as manifested in (|9.2ip . It follows that if $(-;0) is known at a 
fixed 0o then it is known for all = 0o + 2'kNoj + 2tcM with M G 7L d . Then if in addition 
(l,k>) is nonresonant it follows by continuity that it is known for all on the torus. 

Theorem 9.2 Consider a fixed torus Jo off orbital resonance and assume that a O exists 
such that S(0o) 7^ 0- TTien /or euen/ real A and a// G M d , a($(-; 0), A) exists and is 
continuous in 0. Moreover, for all <fi, A($(-; 0)) = A($(-; 0o))- 

Proof: Since $(• + 27rA^; o ) is quasiperiodic, it is easy to see from Definition 4.1 that 

a($(- + 2ixN; o ), A) = exp(iA27rA0a($(-; 0o), A) (9.22) 
and the basic identity ()9.2H) gives 

a($(-; 0o + 2vrA^), A) = exp(iA27riV)a($(-; O ), A)$ T (2vrA^; O ) . (9.23) 

Therefore 

A($(-; O )) = A($(-; O + 2itNu + 2vrM)) , (9.24) 

for all N G Z, M G Z d , where we also used the fact that $(#; O ) is 27r-periodic in O . Thus 
the spectrum of the principal solution matrix on the set 

D' : = {0o + 2iiNuj + 2nM : N eZ,M eZ d } (9.25) 

is the same as the spectrum of the principal solution matrix at 0o- Because (l,u;) is non- 
resonant, D' is dense in M. d . Now fix A and let h(<fr) := a($(-; 0), A) for G D' and assume 
h(<p ) ^ so that A G A($(-; 0o))- Then 

< |/i(0 o )| = |/i(0 o + 2irNu)^(27rN; O )| < v^Ofo + 27riVu;)| 

= V3|/i(0 o + 27rA^o; + 27rM)| , (9.26) 

where we used fl9.23j) for the first equality and where for the second inequality we used the 
fact that 

\X<S>(2ttN; o )| < V3\X\ , (9.27) 

which follows from the SO(3) nature of $. As always, | • | denotes the Euclidean norm, i.e. 
\X\ := y^XnX^ + Xi 2 Xl 2 + ••• + X33X33. If h is continuous and 27r-periodic and defined 
all over M. then, since D' is dense in R d } \h(<j>)\ > for all and thus A($(-; O )) C A($(-; 0)) 
for all 0. Conversely if, for a given 0, |/i(0)| > then, by exchanging the roles of and 0o, we 
obtain |/i(0 o )| > so that A($(-;0 O )) D A($(-;0)) for all 0. Thus A($(-;0 O )) = A($(-;0)) 
for all and it remains to show that h is defined and continuous on R d . The 27r-periodicity 
of h then follows immediately. We first state the following lemma. 
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Lemma 9.3 For fixed A and nonresonant (l,u;), let 

$ g>A (0) := - ! ®(6;(j))exp(-i\6)d6 (9.28) 
Q Jo 

converge uniformly for all in the set D' of ()9.25j) as the positive integer q — > oo. Then <3> 9]A 
converges uniformly on M. d . In particular a($(-; 0), A) exists for all in M d and is continuous 
in 0. 

Proof of Lemma 9.3: The space Y of bounded (w.r.t. Euclidean norm) functions g : M. d — > C 9 
is a complex normed space w.r.t. the norm supj, gK d \g(4>)\ and obviously $ 9)A is a sequence 
in Y . Because C 9 is complete, so is Y (see for example jlHl Section 7.1]). Thus to show that 
$ 9; a converges uniformly on IR d it suffices to show that it is a Cauchy sequence in Y, i.e. for 
all positive 5 there is a positive integer m such that for all integers with j,k > m we have 

sup <$■ (9.29) 

We also observe that 

sup |%, A (0) - $m(0)| = sup \$j,x{<f>) - ®kA<f>)\ > (9-30) 

which follows from the continuity of \<&j t \(<J)) — &k,\{4>) \ in and from D' being dense in 
Equation ()9.30|) implies that ()9.29|) is equivalent to the statement 

sup|% A (0)-$ fejA (0)|<5. (9.31) 

Clearly ()9.31|) holds because, by assumption, $ giA (0) converges uniformly on D' . Thus $ 9i a 
converges uniformly on M d , so that a($(-; 0), A) exists for all in M d and is continuous in 0. 

□ 

To complete the proof of Theorem 9.2 we now show that the conditions of Lemma 9.3 
are fulfilled, i.e. that <& 9)A converges uniformly on D', for every A. Clearly we have 

a($(-; 0o + 2nNu), A) = lim $g, A (0o + 2nNu) 

= exp(iA27riV)a($(-; O ), A)$ T (27riV; O ) (9.32) 

and to show that the limit in ()9.32|) is uniform on D' we estimate 

|$ g , A (0o + 2itNuj + 2nM) - exp(zA27riV)a($(-; O ), A)$ T (27riV; O )| 
= |$ 9i a(0 o + 2ttNlu) - exp(zA27riV)a($(-; O ), A)$ T (27riV; O )| 

i p2nN+q 

= | - / $(0; o )$ T (27riV; O ) exp(-zA(fl - 2-nN))dB 

Q J2ttN 

- exp(iA2vriV)a($(-; O ), A)$ T (2vriV; O )| 

i i>2nN+q 

< \/3|- / $(0;0 o )exp(-a0)ci0-a($(-;0o),A)| . (9.33) 

9 </27r7V 

51 



By using the first equality of (|9.32|) with N = and by noting that $(•; 0o) is bounded we 
have, for all N, 

i f-2-nN+q 

lim |- / $(0; O ) exp(-iX9)de - a($(-; O ), A) | = . (9.34) 

Moreover because $(•; </> ) exp(— zA •) is almost periodic it follows (see jUJ Chapter 3]) that 
the convergence in (|9.34|) is uniform on the domain Z of N. Hence ()9.33|) implies that Q q \ 
converges uniformly on D' . □ 

Thus we have proved, under the conditions of Theorem 9.2, that A(<3>(-; 0o)) is indepen- 
dent of 0o- 
Remark: 

(3) Remark 2 shows that under the conditions of Theorem 9.1 the spin tune can indeed 
be discovered from a spectral analysis of the spin flow for arbitrary 0o- in particular, 
since the spin motion is quasiperiodic and the torus is well tuned, the spectrum has 
at most countably many elements and A($(-;0o)) will consist of sharp "lines" which 
can then be "measured". Moreover, off orbital resonance and under the conditions 
of Theorem 9.2, all A($(-;0 O )) are equal. However, if the torus is ill-tuned but the 
spin motion is quasiperiodic (so that the spectrum A($(-;0 O )) will consist of sharp 
lines), we expect that the union of the spectra A($(-; 0o)) over the torus will contain 
uncountably many elements. The models in Remarks 13 and 14 in Section 6 provide 
examples. Note that in the absence of quasiperiodicity or even almost-periodicity, as 
for example in the Moser-Siegel model of Section 8, there may be difficulties in even 
computing the spectrum. In practice, the spectrum can be obtained by tracking three 
mutually orthogonal spins along an orbit and storing their values at each of a very 
large number of turns before applying a discrete Fourier transform to the data (a well- 
known way of doing numerical Fourier analysis can also be found in |53j). But since 
this spectrum can be very dense, it can be difficult to identify the spin tunes. Thus it 
would be useful to begin with small amplitudes, i.e. tori with small J, and look for a 
uj s close to that for the closed orbit, which can be calculated as in Section 3. The spin 
tune at higher amplitudes could then be identified by continuing away from the closed 
orbit. We return to this theme in Section 10. Finally, tracking a spin trajectory which 
is parallel to an ISF would give a spectrum without u s . 

(4) Of course, there might be cases where the dependence of A($(-; 0o)) on O leading to 
the ill-tuning mentioned in Remark 3, is very weak so that the sharp lines are just 
broadened. Intuition suggests that in such cases, spin-orbit resonance-like phenemena 
might still be expected. But, of course, the definitive conditions under which such 
phenomena could occur, will only become clear by careful analysis. 

10 Discussion and Conclusion 

In the foregoing sections we have presented a thorough step by step account of the circum- 
stances under which spin motion may be quasiperiodic on integrable particle orbits and have 
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thereby put previous studies of the concept of spin tune onto a rigorous basis. In particular, 
we considered integrable orbits in Q(l, u>i, u>d] d + 1) and by introducing certain conditions 
(e.g. Diophantine conditions) and assuming the existence of an ISF we obtained conditions 
under which the spin motion is in Q(l, uj\, uia, uj s ; d + 2) where u s is a spin tune. We 
have also shown how, by introducing UPFs, the spin motion can be represented in terms of 
generalized Floquet forms. 

The scenarios covered by our treatment and the relationships between them are summa- 
rized by the Venn diagram of Figure 1. 



The meanings of the domains in Figure 1 are as follows: 

• Inside the black circle: all tori, i.e. for arbitrary real skew-symmetric 3x3 matrices 
A(-,-,Jq) and arbitrary orbital tunes cu(Jo), which are smooth and 27r-periodic in 9 
and 0. 

• Inside the red ellipse: tori which have an ISF 

• Inside the blue ellipse: tori, which at every 0o have a proper UPF (note that for those 
tori every spin trajectory is quasiperiodic) 

• Inside the green ellipse: tori which are well-tuned (see Section 5) 

• Inside the yellow ellipse: tori which are off orbital resonance 

• Inside the pink ellipse: tori with a uniform IFF 




Figure 1: The logical connections between the various scenarios. 
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The numbered circles label specific examples, namely: 

• Example 1: the single resonance model off orbital resonance (see Section 7) 

• Example 2: the torus J defined in Remark f3 in Section 6, for the case of irrational 

uj(J) 

• Example 3: the single resonance model on orbital resonance (see Section 7) 

• Example 4: the torus defined in Remark 13 in Section 6, for the case u — 1 

• Example 5: a torus defined in Remark 14 in Section 6 

• Example 6: the Moser-Siegel model (for certain choices of the parameters u,ai,(X2) - 
see Section 8 

• Example 7: see below 

Analytical solutions for spin motion or the ISF for the ^4s arising in real storage rings 
cannot be obtained and it is not even known whether an ISF exists in general. Nevertheless 
it seems that it usually exists for storage rings of interest. This is supported by a large 
amount of numerical work in which, for the cases studied, it was possible to construct at 
least a very good approximation to an ISF [HI IE1 El EEU CE21 EH EH ESj - Note that in these 
simulations hard edge and some thin lens representations of fields were used so that the A s 
were not smooth. See Remark 16 in Section 6. 

Those studies also included calculations of the spin tune using SPRINT, either beginning 
with the pseudo-w 1 , w 2 -axes or using the SODOM2 algorithm whereby the spin tune is ob- 
tained by solving an eigenproblem for Fourier coefficients in an SU(2) formalism. Numerical 
simulations with both methods suggest that a spin tune normally exists if the torus is off 
orbital resonance. This is the situation considered by Theorem 6.5d which implies that if we 
have a ISF S and if a 27r-periodic unit vector exists which is nowhere parallel to S, then we 
almost always have a spin tune (see Remark 5 in Section 6). For this and the other results 
in this paper it has been convenient to prescribe that the .As are smooth. But the demon- 
stration of ISFs and spin tunes in El M EEU El EH 1221 122] indicates that the smoothness 
condition can often be relaxed. For convenience, in the remainder of this section we will use 
the term ISF in this spirit. We hope in the future to be able to present a treatment of the 
ISF and the spin tune in which the requirement of smoothness is relaxed. 

An example of a situation where an ISF might not exist even in crude approximation, 
namely near 803.5 GeV in the HERA ring, is given in ^3j- If the spin motion is nonquasiperi- 
odic, this case would correspond to example 7 in Figure 1. 

The spin tune is a crucial quantity for characterizing the stability of spin motion. Spins 
behave somewhat like driven oscillators where the driver is the magnetic and electric fields 
along particle orbits. Near spin-orbit resonance, there is potential for marked qualitative 
changes in the spin motion which may then be quite erratic. The special nature of spin-orbit 
resonance is already clear from the fact that the ISF need not be unique at resonance. See 
Theorem 6.4. Moreover, it is clear that it would make little sense to define a spin-orbit 
resonance condition in terms of a UPR v s depending on <ft Q : such a u s would in general 
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take different values for different particles on a torus so that it would be impossible for the 
particles on a torus to be simultaneously on resonance and there would be no enhancement 
of our ability to systematize spin motion. This was the reason for insisting that a torus be 
well-tuned before considering spin-orbit resonance. 

As mentioned in the Introduction and in Remark 5 of Section 5, it is clear that uj s 
can vary with the beam energy and with J. This is confirmed by simulations in SPRINT 
[HI El EE HH EE21 EE I2H E3 , hence the name "amplitude dependent spin tune" . To study the 
dependence of Sj on J (recall Remark 5 of Section 5), it is necessary to choose a "preferred" 
member of Hj. For the calculations with SPRINT the choice is made as follows. 

The spin tune uq G [0, 1) and the corresponding UPF are found on the closed orbit using 
the method outlined in Section 3. Normally z/ is set significantly different from zero in 
order to ensure that the direction of n is close to the "design" direction in spite of the 
presence of the usual misalignments of the ring (7JEEIHI1- Then as in Remark 3 in Section 
9, the preferred spin tune at nonzero amplitudes is selected by requiring that it and the 
corresponding UPF vary continuously with amplitude and reduce continuously to the spin 
tune and UPF on the closed orbit. 

With this procedure in place, simulations in which the ISF and the spin tune are cal- 
culated over a range of fixed amplitudes or energies indeed show that the spin motion can 
become erratic near a spin-orbit resonance. Moreover, in such cases there is a tendency for 
the ISF to become very sensitive to the parameter being varied. With this prescription, one 
also finds that the strongest variations of the ISF occur near low order resonances and that 
high order resonance effects are usually unimportant. For the details of these calculations 
and results, the reader is referred to [HJ El EE EEH EE2 Ed 1221 • These findings are consis- 
tent with perturbation theoretical calculations of the ISF as in JU] 2 . The nonperturbative 
calculations also show that near spin-orbit resonance the spin tune tends to avoid exact ful- 
fillment of the resonance condition. As we have already indicated in Remark 3 of Section 7, 
the spin tune for the single resonance model (see (|7.11Jl ) also avoids the spin-orbit resonance 
condition as o\ is varied through the condition o\ = u. 

In a ring without so-called Siberian Snakes (see below) the closed orbit spin tune u 
usually varies with the beam energy R is then sometimes implied in the literature 

that in a beam with a large energy spread, the synchrotron motion causes the particles to 
oscillate to and fro across spin-orbit resonances as the spin tune u s oscillates. This crossing 
of resonances is then supposed to be the source of the low beam polarization that would 
be seen. However, we have seen that a spin tune uj s on a torus (Ji, J 2 , ^3) is a constant. 
As usual we assign J 3 to synchrotron motion. Thus uo s does not oscillate and there is no 
resonance crossing. Nevertheless, we do expect that a large energy spread can lead to small 
beam polarization. This is explained by the fact that for particles of large enough J3, n 
varies strongly with the ^3. Then the maximum permissible equilibrium beam polarization 
can indeed be small 0. 



Studies of the effect on polarization of real resonance crossing can be found in ^TJ ^1 
El 123 ■ The consequences for the polarization of crossing first order resonances are usually 
quantified using the Froissart-Stora formula j^lj. But in ^T] EH E3 123] ^ ^ s snown that the 

2 Note that in JOj the perturbation expansions calculate n — n . Then the lo s for the closed orbit appears 
in the resonance factors, not the u> s for nonzero orbital amplitudes. 
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Froissart-Stora approach can be generalized to describe the change of the polarization at the 
crossing of higher order resonances too. This is a further illustration of the value of using a 
wisely defined spin tune for identifying resonances and for understanding their properties 

The spin tune can be obtained in SPRINT via the pseudo-u 1 , u 2 -axes or by using the 
SODOM2 algorithm. However, as we have seen in Section 9 the spin tune might also be 
obtained ("measured") by a spectral analysis of the spin motion. This offers an attractive 
alternative. As in the case of the other two methods the preferred spin tune would have to be 
identified among the many spectral lines by matching onto the spin tune of the closed orbit. 
We have also seen in Section 9 that the ISF might be obtainable, if it exists, by spectral 
analysis of spin motion. Spectral analysis may thus lead to a criterion for deciding whether 
an invariant spin field exists. Other criteria for the existence of the ISF are already available 
for stroboscopic averaging and the SODOM2 algorithm. 

We now complete this discussion by mentioning other quantities that have been used in 
attempts to define a spin precession frequencies. 

As explained in Remark 3 in Section 3, for motion on the closed orbit the spin tune can be 
obtained trivially from the complex eigenvalues of the 1-turn principal solution matrix. This 
does not normally work off the closed orbit since $(27r; <po) = U(2iv) exp (J'27ri , s (U)) U T (0) 
and for this the exp(±2-7rw s (£7)) are generally not the complex eigenvalues. So the reader 
will agree that on synchro-betatron orbits the spin tune usually cannot be obtained from 
the complex eigenvalues of the 1-turn principal solution matrix. In fact in general the real 
eigenvector of the general 1-turn principal solution matrix $(# + 2tt; <f)o)§ T (6; 0o) starting 
at 6 is not even a spin trajectory and is not parallel to an ISF [55, p. 27, for example]: the 
term "spin closed orbit" which is sometimes used for the ISF, is inappropriate. Moreover 
in the simulations with SPRINT the sensitivity of the ISF to variation of parameters shows 
no correlation with the spin precession rate extracted from the complex eigenvalues of the 
1-turn principal solution matrix. Of course eigenvalues of 1-turn principal solution matrices 
are easy to calculate but for J ^ they usually have no useful function. 

The model described in |3*2"1 equation 21] provides an example. It concerns orbits which 
are said in [3T] to be exceptional. It is stated there that exceptional orbits are characterized 
by the feature that the spin tune depends on orbital phase and this is discussed in close 
conjunction with Stern-Gerlach forces. As we point out in the Introduction, S-G forces can 
have no practical relevance for understanding spin resonance. The model in |J2j involves 
the single resonance model (see Section 7) and a single thin lens Siberian Snake. A Siberian 
Snake jSHl GUI °f the kind used here is a magnet system that rotates a spin by the angle 7r 
around an axis in the plane of the ring p. 71]. In the notation of Section 7 the parameters 
in (32] are o\ = lo and a 2 V2J = 3/2. At all 4>o the 1-turn principal solution matrix $(27r; <j) ) 
starting at the snake represents a rotation around the vertical by an angle depending linearly 
on 0o- Then the ISF is vertical at the snake at all 0o o n this torus and for this zero-measure 
range of J, the ISF at the snake is the same as the real unit-length eigenvector of the 1- 
turn principal solution matrix. 3 The orbital tune u is arbitrary. The eigentunes extracted 
from the complex eigenvalues of $(27r;0 o ), vary linearly with <j) . Therefore if a particle is 
followed along an orbit, these eigentunes normally change abruptly between one turn and 
the next. Other 1-turn eigentunes would be obtained if other positions around the ring were 

3 Note that in this discussion we allow the ISF to be discontinuous in 6. 
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chosen for 9 = 0. Nevertheless, it is claimed in E2] that the eigentune for $(27r;0o) is 
27r x spin tune. Obviously this is not a spin tune in the sense of our treatment. We have 
taken account of the fact that an eigentune in the SU(2) formalism used in |3T|l32] is one half 
of a corresponding eigentune in our formalism. Although this $(27r; <fio) and its eigentune 
are discussed in association with S-G forces, the inclusion of S-G forces is not necessary to 
obtain either. In fact by exploiting techniques additional to those used in this paper one can 
show that these parameters simply provide an example of ill-tuning and that this case would 
be entered next to example 6 in the diagram of Figure 1. The inclusion of S-G efffects clouds 
the issue. It is not clear from [HH E2J whether these phase dependent eigentunes serve some 
useful function such as indicating the stability of spin motion. We give another example of 
the use of the eigentune of a 1-turn principal solution matrix below. 

The situation with regard to the utility of eigenvectors and eigentunes is more subtle for 
tori on orbital resonance where u\, u>2, ... are rational. The orbit and A are then periodic over 
an appropriate number, M, of whole turns. Thus in analogy with the method in Remark 6 in 
Section 3, there is a possibility of calculating n at 9 = as the unit-length real eigenvector 
of the M-turn principal solution matrix and in general it would be a function of 4>q. The 
imaginary part of the exponent of a complex eigenvalue of this principal solution matrix 
would provide the M-turn advance of the phase of spin rotation around n and this could 
be used to obtain the average 1-turn spin phase advance. This would usually depend on O 
and in such cases it could not be used to define a spin tune. But it could be used to define 
a z/y (see (|fi.2fi|) ). 

So, as is usual at orbital resonance, an ISF can exist in general but normally there is 
no spin tune. We note in passing that the M-turn complex eigentunes are 27r-periodic in 9 
and 0, just like A. Moreover, because eigenvalues of matrices are invariant under similarity 
transformations, the eigentunes are invariant when the starting point for the eigenanalysis is 
shifted along the orbit. Examples of an unwise use of the term spin tune can be found in [5%] 
where the dependence of the multi-turn eigentune on 4> is made explicit. See jSHE2] too. 
Again, a "tune" depending on 0o cannot be used for studying spin-orbit resonance. Calling 
such a quantity a spin tune can create confusion. See Remark 4 in Section 9 also. 

If, on orbital resonance with rational tunes, the torus is ill-tuned, then one can expect 
that either there are more or fewer proper UPRs than one would have on a well-tuned torus. 
Thus a spectral analysis along the lines of that in Section 9, applied to every O on the torus, 
could be a useful diagnostic tool to signal these two cases of ill-tuning. Note that in the 
examples in Remarks 13 and 14 in Section 6 there are too many UPRs. Therefore one might 
expect that on orbital resonance with rational tunes an ill-tuned torus had too many proper 
UPRs. 

Although an M-turn eigentune is in general 0o dependent, we might expect that an 
approximation to the spin tune u s on a well-tuned torus off orbital resonance could be 
obtained by setting the orbital tunes to rational values near to the actual tunes but such 
that M were very large. Indeed this is the essence of a popular perception. In effect, although 
usually not clearly stated, the underlying hope is that for the smooth guide fields of real 
rings and for large enough M, the average 1-turn spin phase advance is only very weakly 
dependent on O and that therefore a good approximation to u s can be obtained. That 
would be consistent with the heuristic expectation that for large enough M the influence 
of the initial O become diluted. Such behavior might also be expected if the high order 
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Fourier coefficients in the cy(6 l , — uo9) in (j6.21J) are very small owing to the smoothness of 
the fields. See |24[ p. 66] for a hint of how the Fourier coefficients come in here. We hope in 
the future to be able to present a rigorous treatment of this approximation. This approach 
has been adopted in [SB] to indicate that off orbital resonance the spin tune is a half integer 
at most J in the model to be discussed in the next paragraph. However, any attempt to 
find an approximate value for a spin tune by using M-turn maps for nearby rational tunes 
should at least be checked for convergence and consistency. See Remark 4 in Section 9 too. 

Further discussion around the topic of rational tunes can be found in jSI] in which the 
nature of the so-called "snake resonances" is studied. These refer to a large loss of beam 
polarization during acceleration in a model in which the spin motion in most of the ring 
is approximated by the single resonance model and the spin motion is stabilized by pairs 
of idealized, i.e. thin lens, Siberian Snakes The Siberian Snakes have the effect of 

fixing u at 1/2 independently of the beam energy jS^l p. 70]. Snake "resonances" occur at 
rational orbital tunes uj for which, in the notation of (|5.1U|) . 1/2 = m + m% uj with odd 
m 1 . In |HJ] it is made evident that at these tunes and for most J ^ there is no amplitude 
dependent spin tune so that one is not dealing with spin-orbit resonances. Moreover at most 
orbital amplitudes the ISF defined there (i.e. without insistence on smoothness) is irreducibly 
discontinuous at some orbital phases [HI]. Of course, since rational tunes correspond to 
orbital resonance, it should come as no surprise, given the content of our paper, that there is 
no amplitude dependent spin tune in this case. We have explicitly chosen our nomenclature 
to be consistent with earlier usage [IB] and thereby contribute clarity to the classification 
of phenomena. It could be that the pathological behavior in this model, namely the large 
loss of polarization, is due to the use of the simplified but singular representation for the 
snake fields. In this connection we note with interest that according to simulations for 
RHIC, the loss of polarization during acceleration is less severe when the simulations are 
carried out with the magnetic fields of real snakes than with the singular fields of thin lens 
snakes jHO] ■ This indicates that predictions from simplified, mathematically singular models 
should be treated with some caution. For recent experimental work around vertical orbital 
tunes corresponding to "snake resonance", see [HI]. In any case the use of the term "snake 
resonance" is a good illustration of the confusion that arises from an imprecise use of the 
concept of spin-orbit resonance. 

A model involving the single resonance model and two thin snakes provides the second 
example of the use of the eigentune from the 1-turn principal solution matrix. Simulations 
described in |59| p. 98-100] show some loss of polarization for all nonzero 2<j\ J during accel- 
eration, even away from the rational orbital tunes associated with snake "resonances" . It is 
implied there that this loss stems from the fact that during the acceleration, the eigentune 
from the 1-turn principal solution matrix, which in [2^1 is called the "perturbed spin tune" , 
oscillates to and fro across a spin-orbit resonance as the orbital phase advances from turn to 
turn. But in this example the 1-turn eigenvector of the principal solution matrix is usually 
not even parallel to the ISF. So it is even more difficult to imagine that the "perturbed spin 
tune" can characterize long term spin motion. Thus an alternative picture for the loss of 
polarization must be sought. 

As mentioned in Remarks 1 and 4 in Section 7, for the straightforward single resonance 
model a spin tune does exist on orbital resonance. 

This completes our discussion of notions of spin precession frequency. 
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We now conclude by summarizing the main message. This is, that, by just employing the 
Lorentz force and the T-BMT equation and with the help of the concept of quasiperiodicity, 
we are able to provide a rigorous and broad treatment of the concepts of proper uniform 
precession rate, spin tune and the invariant spin field by common methods used for ordinary 
differential equations. This allows us to focus on the main phenomena without the distraction 
of perturbations such as noise, collective effects and synchrotron radiation. In principle 
they can be included by using perturbation theory. We have discussed putative Stern- 
Gerlach forces in the Introduction. The advantages of a clear, universally accepted and useful 
definition of spin tune has been made evident by the examples in the foregoing paragraphs. 
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